Extending Luttinger's theorem to Zi fractionalized phases of matter 
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Luttinger's theorem for Fermi liquids equates the volume enclosed by the Fermi surface in momen- 
tum space to the electron filling, independent of the strength and nature of interactions. Motivated 
by recent momentum balance arguments that establish this result in a non-perturbative fashion [M. 
Oshikawa, Phys. Rev. Lett. 84, 3370 (2000)], we present extensions of this momentum balance 
argument to exotic systems which exhibit quantum number fractionalization focussing on Z2 frac- 
tionalized insulators, superfluids and Fermi liquids. These lead to nontrivial relations between the 
particle filling and some intrinsic property of these quantum phases, and hence may be regarded as 
natural extensions of Luttinger's theorem. We find that there is an important distinction between 
fractionalized states arising naturally from half filling versus those arising from integer filling. We 
also note how these results can be useful for identifying fractionalized states in numerical experi- 
ments. 
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I. INTRODUCTION 



The last two decades have witnessed the experimen- 
tal discovery of several strongly correlated materials that 
show properties strikingly different from that expected 
from conventional theories based on Landau's Fermi liq- 
uid picture. These include the high temperature cop- 
per oxide superconductors, heavy fermion systems near 
a quantum critical point, and, more recently, the cobalt 
oxide materials. Interesting correlated quantum phases 
are also likely to emerge in the near future from ongoing 
experimental efforts in the area of cold atoms in optical 
lattices. It has then become imperative to theoretically 
investigate quantum phases of matter that differ funda- 
mentally from the standard paradigm. Indeed, in such a 
search for new theoretical models, it would be useful to 
know if general principles place constraints on the possi- 
ble quantum phases. Here, we will explore in detail the 
consequences of one such constraint, arising from mo- 
mentum balance, which will be applicable to interacting 
many body systems on a lattice. This argument was first 
applied to the case of one-dimensional Luttinger liquids^, 
where it relates the Fermi wavevector kp to the particle 
density. It was later extended to Fermi liquids^ in di- 
mensions D > 2, where it leads to Luttinger's theorem-, 
relating the filling fraction to the volume enclosed within 
the Fermi surface on which the long lived Fermi-liquid 
quasiparticles are defined. Here, we will apply the same 
line of argument to a variety of different phases in spatial 
dimensions D > 1, and the constraint we obtain in this 
way may be viewed as analogues of Luttinger's theorem 
for these phases. In all cases, the filling fraction (number 
of particles per unit cell of the lattice) is fundamentally 
related to some intrinsic property of the phase. 

The momentum balance argument, introduced by Os- 
hikawa for Fermi Liquids-, proceeds as follows. Consider 



a system of interacting fermions at some particular filling 
on a finite lattice at zero temperature. Periodic boundary 
conditions implies that the lattice has a torus geometry; 
imagine introducing a solenoid of fiux in one of the holes 
of the torus, and adiabatically changing its strength from 
zero to 277. The crystal momentum imparted to the sys- 
tem can then be caculated in two different ways. First, 
in a trivial fashion that only depends on the filling and 
is independent of the quantum phase the system reaches 
in the thermodynamic limit, and second, in way that de- 
pends essentially on the quantum phase of the system. 
Consistency requires the equality of the these two quan- 
tities - which leads to the nontrivial conditions on the 
quantum phase. Essentially, each consistent quantum 
phase has its own way of absorbing the filling dependent 
crytal momentum that is generated in this process - as 
mentioned, in the case of the Fermi liquid this leads to 
Luttinger's theorem. 

Here, we begin by applying this argument to the case 
of (bosonic) insulators at half filling, where the system 
in the thermodynamic limit necessarily acquires an en- 
larged unit cell (through broken translational symmetry 
or a spontaneous flux) or develops topological order. For 
the latter case, the momentum balance argument fixes 
the crystal momentum of the degenerate ground states 
in the different topological sectors. A useful side result 
of this analysis will be a general prescription to distin- 
guish between between a Z2 fractionalized insulator (or 
spin liquid) and a more conventional translation sym- 
metry broken state, which is useful when the the order 
parameter for the translation symmetry breaking is not 
obvious. This is relevant for numerical studies on finite 
sized spin systems that search for fractionalized spin liq- 
uid states^. 

Next, we apply the same methods to the case of exotic 
fermi liquids {FL* phases) proposed recently in Ref. |a. 



This is a phase which has conventional electron like exci- 
tations near a Fermi surface, but also posesses topological 
order and gapped fractionalized excitations. The ques- 
tion of interest here is whether the Fermi surface in these 
systems 'violate' Luttinger's theorem (given that these 
phases are not continuously connected to the free elec- 
tron gas this is of course not prohibited by Luttinger's 
proof'^), and if so whether there is a generalization of 
Luttinger's theorem that can accomodate these cases as 
well. Indeed, applying the momentum balance argument 
to the FL* phases we find that while Luttinger's Theo- 
rem is violated by these Fermi volumes, this violation is 
not arbitrary but is constrained to be one of a few pos- 
sibilities which is determined uniquely by the pattern of 
fractionalization. 

Finally, we apply these arguments to the case of neutral 
superfluids. For conventional superfluids we argue that 
this leads to a constraint on the Berry phase acquired 
on adiabatically moving a vortex around a closed loop, 
by relating it to the boson filling. Loosely speaking this 
is the quantity that determined the Magnus 'force' on a 
moving vortex. Since this relation between the Magnus 
force and the boson filling is obtained using the same mo- 
mentum balance argument that leads to Luttinger's theo- 
rem when applied to a Fermi liquid, it may be viewed as a 
'Luttinger' theorem for superfluids. Alternatively, since a 
similar relation between boson density and Magnus force 
is known for superfluids with Galilean invariance^'-, this 
may be viewed as an extension of those results to the 
case of lattice systems. In contrast, while Galilean in- 
variance also constrains the zero temperature value of 
the superfluid stiffness, that constraint does not survive 
the inclusion of the lattice. In order to further bring 
out the similarity of this relation in superluids to the 
Luttinger relation, we consider fractionalized superfluid 
phases SF*, (related to the exotic superconductor SC* 
studied in refiS), or equivalently superfluid analogues of 
the fractionalized Fermi liquid phases. We show that 
they too violate the conventional relation between vor- 
tex Berry phase and boson filling in exactly the same way 
that Luttinger's theorem is violated in FL*. We discuss 
caveats in the relation between the vortex Berry phase 
and the boson filling in conventional and fractionalized 
superfluids which make the above relation less rigorous 
at the present time than the analogous relation for Fermi 
liquids and insulators. 

The relation between vortex Berry phase and boson 
filling in lattice superconductors can lead to surprising 
conclusions in some cases. For example, consider a con- 
ventional superfluid (where the bosons are Cooper pairs 
of electrons) obtained by doping a band insulator ver- 
sus another conventional superfluid obtained on doping 
a proximate Mott insulator. One may imagine that only 
the doped charges participate in the superfluidity - in- 
deed this is roughly what is expected for a quantity like 
the superfluid stiffness (although it is not strictly con- 
strained in these lattice systemsiS). However, a result 
of the discussion below will be that the Berry phase ac- 



quired by a vortex in this system arises from counting all 
charges in the system (and not just the charges doped 
into the Mott insulator). In this sense at zero tempera- 
ture all particles participate in the superfluidity. In con- 
trast, an exotic superfluid phase SF* can display a phase 
where only the doped charges contribute to the Berry 
phase. 

A recurring theme throughout this paper will be the 
distinction between exotic phases obtained from a corre- 
lated 'band' insulator (i.e. one that has an interger filling 
per site in the case of bosons) and which could form a 
translationally invariant conve ntional insulating state; 
versus those obtained by from an exotic phase at half 
filling. For example, if we consider featureless Z2 frac- 
tionalized insulators at integer and half integer filling, 
then at low energies they are be described by 'even' and 
'odd' Z2 gauge theories respectively (in the terminology 
of Ref. 11). The different ground state topological sectors 
of the odd gauge theory can in certain geometries carry a 
finite crystal momentum, while the the crystal momenta 
associated with the different ground state sectors of an 
even gauge theory are always zero. This distinction per- 
sists if these phases are doped to obtain exotic Fermi liq- 
uids and superfluids. The distinction is especially strik- 
ing in the case of FL*^^^ where a Fermi volume that vi- 
olates Luttinger's theorem arises. In contrast, FL^^^^ 
obeys Luttinger's theorem but is nevertheless an exotic 
phase. A similar distinction will apply to the exotic su- 
perfluids - in that case the relation between the Magnus 
'force' on a vortex and the filling is the regular one for 
^-^even t>ut is unconvcntional in the case of SF*^^. 

An essential ingredient in the following arguments will 
be the evolution of a quantum state under fiux insertion. 
While this recalls the argument of Laughlin^^ for the in- 
teger Quantum Hall effect, there is an important distinc- 
tion that must be noted. In the case of Laughlin's argu- 
ment and a similar argument for the forces on superfluid 
vortices given by WexleBi"^, the conclusions are derived by 
keeping track of the change in energy during the process 
of flux threading. More recently, rigorous energy count- 
ing arguments for charge and spin insulators have been 
made by Oshikawcii^ and Hastingai^. In contrast here 
we will follow Ref. 2 and rather keep track of the change 
in crystal momentum during the flux threading process, 
which will allow us to derive a different set of rather gen- 
eral conclusions that apply to a variety of phases. We 
also note that while the subject of Magnus force on a 
vortex at finite temperatures, in the presence of quasi- 
particle or superfluid phonon excitations, has been the 
subject of much lively debate (see for example Refs.fig), 
our arguments will only apply to the case of zero temper- 
ature and hence cannot address any of the issues under 
debate. 

The layout of this paper is as follows. Due to the length 
of the paper we give in Section ^ a simplifled overview 
of all results and a brief indication of the method used. 
Then, we pass to the technical details and in section IhTI 
review the momentum counting procedure which will be 



applied to all the phases. Then, we consider conventional 
insulators in section IIVI and Z^ fractionalized insulators 
in section fusing the momentum balance argument and 
discuss how they may be unambiguously distinguished 
in numerical experiments in section IVII We then re- 
view the momentum balance argument for conventional 
Fermi liquids in section IVlII which leads to Luttinger's 
theorem, and see how this is modified ina systematic way 
when applied to Z^ fractionalized Fermi liquids in section 
IViiil Next we apply these arguments to conventional, 
neutral superfluids and Z2 fractionalized supcrflui ds, in 
Sections Hxl and Ixl respectively. We conclude with some 
observable consequences that arise directly from these 
considerations. 



II. OVERVIEW OF THE MOMENTUM 
BALANCE ARGUMENT 



In this section we summarize the results of the momen- 
tum balance argument applied to different phases. While 
the detailed arguments leading to these results are con- 
tained in the following sections, the results themselves 
are easily stated, which is done below along with some 
heuristic supporting arguments. 

Consider the system in a cylindrical geometry, as 
shown in Fig. ^ with dimensions L^^ Ly (integers), and 
a total of N particles (bosons or spinless fermions). Now 
imagine adiabatically threading a flux of 2n through the 
center of the cylinder - the particles are assumed to cou- 
ple minimally to this flux with a unit charge The total 
momentum imparted to the system can be calculated us- 
ing Faraday's law F^ = —j^-^ and integrating this force 
over time leads to the change in momentum 






(1) 



which is only defined modulo 27r since the system is on 
a lattice. A more rigorous calculation in the following 
sections arrives at the same result. A similar procedure 
performed along the perpendicular direction (for which 
it is convenient to think of the system living on a torus) 
yields the other component of the crystal momentum 



27r 
AP„ = — iV. 



V 



(2) 



This is the result of trivial momentum counting - we now 
consider how this additional crystal momentum is acco- 
modated in the various different phases. For later pur- 
poses, we will define the "filling" v = N / L^Ly, which is 
the number of particles per unit cell. 
Insulators: Consider the evolution of the ground state 
of a many body system under the adiabatic insertion of 
a 27r flux. For an insulator, it may be easily verified 
that the final state must have an energy equal to that 
of the ground state, in the thermodynamic limit. That 
is, the system ends up either in the original state, or 




FIG. 1: Schematic figure showing flux threading in a cylinder 
geometry, with flux $0 = hc/Q. 



in a degenerate ground state. This is just because the 
change in energy under fiux threading is related to the 
average current (J) = ^. Since in the insulating state 
the current must vanish in the thermodynamic limit, the 
energy of the final state must equal that of the initial 
state. 

If the insulator is at integer filling, i.e. N — vL^Ly, 
with the filling v such that i^ € integer, then we have 
(APa; = 2TnyLy, APy = 2ttvLx), and so AP^ = APy = 
O(mod 27r). This is compatible with the system having 
a unique ground state which it returns to at the end of 
the 27r flux threading. This is the conventional feature- 
less insulating state (band insulator for Fermions, integer 
filling Mott insulator for bosons) - although in principle 
more complicated states are possible at integer filling as 
well. 

The case of insulators at non-integer filling is more in- 
teresting. For definiteness, consider bosons at half filling 
(1/ = 1/2). In order that the total number of particles 
be an integer, we need the total number of sites L^ x Ly 
to be even. If we first consider the case Ly odd and L^ 
even, under flux insertion in the geometry of Fig. ^ we 
will have AP^; = 7r(mod27r). Therefore, the initial and 
final states (which we have argued to be denenerate in 
the thermodynamic limit) must differ in crystal momen- 
tum and hence one is forced to conclude that the ground 
state is at least doubly degenerate in this even x odd ge- 
ometry. Such a degeneracy can result from one of two 
different reasons (we assume that time reversal symme- 
try is not spontaneously broken and the special case of 
TT-fiux is discussed in Section ^3 D). First, the system 
may be heading towards translation symmetry break- 
ing in the thermodynamic limit. In this case we can 
form the symmetric and antisymmetric combinations of 
the two ground states which clearly transform into each 



other under a unit horizontal translation. Translation 
symmetry breaking implies that there is a local operator 
that can distinguish between these two states. For ex- 
ample, if the system is heading towards a CDW state, 
(eg. with a stripe pattern with the charge on alternating 
columns in Fig. QJ, then the relevant local operator is 
simply the charge density, which would distinguish these 
two states as being translated versions of one another. 
It may of course happen that the relevant local operator 
is less obvious (eg. bond centered charge density) but 
nevertheless in principle this distinction between the two 
states can be made with some local operator. It may 
however happen that there exists no local operator that 
can distinguish these two states. Then, the system will 
appear perfectly translation symmetric, although it is an 
insulator at half filling. Indeed this is precisely the prop- 
erty of the RVB spin liquid state proposed for spin 1/2 
lattice systems with one spin per unit cell - which can 
also be cast in the language of the above discussion on 
identifying the spins with hard core bosons. Therefore, 
these degenerate ground states can only be distinguished 
via a non-local operator. This is called topological or- 
der, where degeneracies arise on spaces with non trivial 
topology that are related to creating a topological excita- 
tion which is highly nonlocal in the original variables. In 
the following sections we consider one concrete theoreti- 
cal realization of this scenario - the case of a deconfined 
Z2 gauge theory coupled to bosons carrying one half the 
elementary unit of charge. The translationally symmet- 
ric state at half filling may be roughly visualised as a 
uniform state with a half charge at each lattice site. The 
degenerate ground states correspond to the topological 
degeneracy of the theory on a cylinder - which is relate 
to the presence or absence of an Ising flux (vison) in the 
hole of a cylinder. In the following sections we explic- 
itly demonstrate that threading a 2tt flux induces such 
a topological excitation and causes the ground state to 
evolve into this distinct topological sector. Thus, the way 
these topologically ordered states accomodate the crys- 
tal momentum imparted to the system on flux threading, 
despite being translationally symmetric, is by creating a 
vison excitation in the hole of the cylinder which then 
carries the appropriate crystal momentum. Later we will 
draw a distinction between Ising gauge theories where 
the vison carries a crystal momentum (odd gauge theo- 
ries) and those where it does not carry momentum (even 
gauge theories). 

Fermi Liquids: The original application of the momen- 
tum balance argument was to the case of Fermi liquids 
in Ref. [^ This argument is reviewed in the following 
sections - here we just note the main points. If we begin 
with spinless electrons at a filling v, then the flux thread- 
ing excites quasiparticles around the Fermi surface. The 
total crystal momentum carried by these excitations can 
be converted into an integral over the volume enclosed by 
the Fermi surface, which leads to the relation between the 
filling (which enters the trivial momentum counting) and 
the Fermi volume. For an appropriately chosen system 



size, these lead to the relation: 

Vfs 



{2ny 



(3) 



where Vps is the Fermi volume, and p is an arbitrary 
integer which represents the filled bands. This is just 
Luttinger's Theorem - in Ref. Q it was also applied to 
Kondo lattice models where it yields the large Fermi sur- 
face expected in the Kondo screened phase. 

One can now ask the reverse question — given a phase 
whose low energy exciations are electron like Landau 
quasiparticles, does this phase necessarily also satisfy 
Luttinger's theorem? From the above momentum bal- 
ance argument it is clear that in order to violate the 
Luttinger relation there must exist an alternate sink for 
the momentum. From our previous discussion of topo- 
logically ordered states, it is clear that if topological or- 
der coexists with Fermi liquid like excitations, then the 
momentum balance can be satisfied with a non-Luttinger 
Fermi volume. In fact the specific case of an exotic Fermi 
liquid with Z2 topological order (FL* phase) was pro- 
posed in Ref. 6 in the context of the heavy fermion sys- 
tems. This phase has low energy excitations identical to 
that of a Landau Fermi liquid of electrons, but also a 
gapped Ising vortex excitation. The Luttinger relation 
relating the filling to the volume enclosed by the Fermi 
surface of these quasiparticles can then be violated, but 
in a very specific way. Flux threading creates a Ising vor- 
tex which carries crystal momentum n (in an odd x even 
system), while the remaining momentum is absorbed by 
the quasiparticle excitations. This leads to the modified 
Luttinger relation: 



Vfs 
(2^)2 
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(4) 



where p is an arbitrary integer representing filled bands. 
Note the crucial difference from Luttinger's relation in 
Eqn|31 that arises from the extra factor of i. Clearly, 
this is related to the fact that a translationally invari- 
ant insulator is possible at half filling, where the Fermi 
volume can shrink to zero. Thus, the difference from 
the original Luttinger relation is very specific, i.e. the 
filling of exactly half a band, for the case of Z2 fraction- 
alization. We have commented earlier on the difference 
between 'odd' vs. 'even' Z2 gauge theories, again this 
distinction is crucial here and moreover is not directly 
set by the filling as it was for the case of the translation- 
ally symmetric insulating states. The above violation of 
the Luttinger relation only occurs in the case of the 'odd' 
gauge theory, FL*^^^. 

Superfluids: Finally, we consider the case of neutral 
superfluids. Here threading a 27r flux clearly inserts a 
vortex through the hole of the cylinder. This can also be 
visualized as creating a vortex in the superfiuid at the 
bottom of the cylinder and dragging it all the way to 
the top. Clearly such a vortex will experience a 'Mag- 
nus force' in the direction perpendicular to its motion. 



Let us ignore for a moment the lattice and calculate the 
momentum imparted by this force Fm = liraMV x i, 
where v is the velocity and aM a constant that fixes the 
Magnus force. The total momentum transferred to the 
system is then independent of the details of the vortex 
motion and depends only on its net displacement — this 
yields AP^: = 2T:aMLy. Equating this to the momen- 
tum obtained from trivial momentum counting JQ and 
reintroducing the lattice heuristically by allowing the the 
momentum to change in arbitrary integer multiples of 27r 
we have: 



ly = aM +P 



(5) 



where p is an arbitrary integer. Thus, the fractional part 
of Um is completely determined by the boson filling i/. 
This can be viewed as the analogue of Luttinger's the- 
orem for bosonic system, since it is obtained using the 
same line of argument. It can also be viewed as an exten- 
sion of the well known equivalent result for Galilean in- 
variant superfluidsi (where the Magnus coefficient is the 
boson density) to the case of lattice superfluids. While we 
have been interpreting aM above in terms of the 'Mag- 
nus force' clearly this is not a well defined concept in a 
lattice system. In fact, the property that is sharply fixed 
by aM is the Berry phase acquired by a vortex on adia- 
batically taking it around a big loop of size A/" plaquettes, 
which will be shown to be 2-KaM-N' ■ This is related to the 
well known relation^ in Galilean superfluids between the 
Magnus force and the Berry phase acquired by a vortex. 
Again, one can ask if the relation is Eqn.© can be 
violated in any kind of superfluid. Indeed, topologically 
ordered superfluid states SF* can be defined in complete 
analogy with FL* . For the particular case of an exotic 
superfiuid state SF* with Z2 topological order'^, there 
exists in addition to the usual vortex excitation, an Ising 
vortex excitation as well. Threading 27r flux then creates 
both a regular vortex and an Ising vortex - the latter 
can carry crystal momentum tt (in the phase SF*^^^^ - in 
which case the remaining momentum is associated with 
the 'Magnus force' on the superfluid vortex. Then, in 
this case as well, the Magnus coefficient aM, associated 
with the vortex Berry phase satisfies: 



'^ = 2 + "M + P 



(6) 



where p is an arbitrary integer. Thus, the 'Luttinger 
relation' for a conventional superfiuid in Eqn.ljSJ is vio- 
lated, in exactly the same way that FL*^^^^ violates the 
Luttinger relation for conventional Fermi liquids. 



III. TRIVIAL MOMENTUM COUNTING 

Consider a system of N particles, each with charge 
Q, living on an L^ x Ly lattice wrapped into the form 



of a torus with periodic boundary conditions along both 
directions. The main result of this section is that if one 
adiabatically threads flux hc/Q through one of the holes 
of the torus, the crystal momentum difference between 
the initial and final state is: 

Pf - P,^ 2TTN/L^{-mod 2tt) = 2TnyLy{-mod 27r) (7) 

where v = N /{LxLy) is the charge density in units of Q 
(the "filling"). This result is independent of the even- 
tual quantum phase of the system in the thermodynamic 
limit, and we will refer to it as the "trivial" counting. 
Although this has been shown in Ref. |2, we include a 
derivation here for the sake of completeness, and to fix 
notation. 

The Hamiltonian for an interacting set of particles 
(fermions or bosons) coupled to an external vector po- 
tential can be written down as Ha = Ka + V^[?^]j with 
the kinetic energy Ka, in the presence of a vector poten- 
tial Aij , given by 



Ka 



2^ 



/ .e-'Q'4,,AcotH ^^^^^ 



(8) 



Bj creates a particle carrying charge-Q at site i. The 
interaction term V[n\ depends only on the density of the 
particles. Both i^ and V[n] are invariant under lattice 
translations. For simplicity of presentation, we will as- 
sume that Uj only connects nearest neighbor sites on a 
square lattice, with unit lattice spacing. 

To thread a unit flux $0 = hc/Q through the hole of 
the torus, say along the —y axis as in Fig. ^ we can 
choose a uniform gauge in which Ai^i^x — ^^{t)/Lx and 
Aij = for other links, and adiabatically increase <I>(i) : 
-^ hc/Q. The state we reach on flux insertion can of 
course be written as \'^{T)) = Ux\'^{0)) with the unitary 

time evolution operator Ut — "^ cxp I — i J^, HA{t)dt] 

where Tt is the time-ordering operator. In the final 
Hamiltonian, the vector potential corresponds to flux 

$(r) = $0- 

Clearly, the initial and flnal wavefunctions, as well 
as the Hamiltonian, transform under gauge transforma- 
tions. Thus, since the final Hamiltonian includes a unit 
flux quantum, we need to flx a gauge in order to con- 
sistently define the crystal momentum of a state as the 
eigenvalue of the unit lattice translation operator acting 
on the state and to compare it for the two states. We 
pick a gauge such that Aij = in the initial as well as 
the final Hamiltonian — this Hamiltonian is denoted as 
Hq. In this case, for a threaded fiux $0, we need to make 
a unitary gauge transformation 



Ha{T) ^ UgHa{T)U^^ ^ Ho 



with the operator 



Kg = exp(i— ^ Xifii). 



(9) 



(10) 



The final wavefunction in this gauge is, in obvious no- 
tation, l^y) — UoUTY^i)- To compute the crystal mo- 
mentum of this state^ we must act on it with the unit 
translation operator T. This defines the initial and final 
crystal momenta, Pi,Pf through T\'^i) = cxp{—iPi)\'^i) 
and T|*/) = exp{-iPf)\'^f). Translating the final state 
we find 

= {fUGf-^)UTe-'''^\^^), (11) 

since the operator IAt commutes with the T as the time- 
dependent Hamiltonian is translationally invariant in the 
uniform gauge. At the same time, it is straightforward 
to show that 



where the average current in units of Q is 



TUgT-^ = exp(-i27riV/L^)Z^G. 



(12) 



It is then clear that Pf = Pi + 2TrN/Lx{raod 27r), or 
defining the filling v — N/{LxLy), the change in crystal 
momentum is Pf — Pi = AP — 2T:i>Ly{'niod 27r). 

It is essential for this argument to go through that one 
has a conserved U{1) charge, this permits us to couple the 
charge to an inserted solenoidal flux. One can easily gen- 
eralize the argument to cases where the charged particles 
carry spin and are coupled to spins fixed to the lattice 
such as in a Kondo lattice modelSi. In this case, one can 
thread a flux which couples to a single component of the 
spin of the charged carriers, and eliminate the vector po- 
tential using a unitary transformation which acts on the 
charged particles as well as the fixed spins. 

As mentioned earlier, the result above has been de- 
rived without any assumption about the thermodynamic 
phase of the system. Such an assumption is important 
for counting the momentum in a second independent 
way, which provides constraints on the various quantum 
phases of the system and we turn to this in the remain- 
ing Sections. For convenience of notation, we will set 
h^c=l in most places. 



IV. CONVENTIONAL INSULATORS 

A. No Broken Symmetry 

Consider a conventional insulator with a unique 
ground state and a nonzero gap to current carrying exci- 
tations. Under adiabatic flux threading, since the Hamil- 
tonian is time-dependent, the rate of change of energy is 
given by (dH/dt) = -J2^j {Jijit))dA,j{t)/dt where the 



current operator J.;. 



iQUjiBJB^e-^^'^'^^^ - h.c.) for 



the Hamiltonian with kinetic energy as in|Sl Let us as- 
sume a linear rate of change of Aij (t) (for j = i+x) over a 
time interval T for threading one flux quantum <I>o = 27r, 
i.e. the electric field Eij = —dAij/dt = {2tt/QLxT)x is a 
constant over the interval T. The total change in energy 
is thus 

i-T 



/ = l/(?igL,r) V / dt{J,.,+i{t)). (14) 



Clearly / = in an insulator in the thermodynamic 
limiliii — there is no current flow, and thus 5E = 0! 

However, if we thread one flux quantum into the sys- 
tem it can be eliminated using a gauge transformation 
which leaves the spectrum invariant, as is well-known 
and was shown in the previous Section. Since the sys- 
tem has a unique ground state with a charge gap, and 
5E = 0, this means the final state and the initial state 
in the Aij = gauge must be the same. Clearly, there 
is no change in crystal momentum on threading flux $0, 
which implies 



2TTi/Ly = O(mod 2tt) 



(15) 



dE = I dt{dH/dt) = 27rhl, 
/o 



(13) 



for any Ly. This is only possibly if v is an integer. Thus 
we arrive at the result: a conventional insulator with 
a unique ground state (i.e., no broken symmetry) and 
a nonzero gap to charged excitations is only possible at 
integer fillingiiSiiS 



B. Conventional insulator with broken 
translational symmetry 

Imagine tuning the interaction V[n\ in the above 
Hamiltonian in Eq. |S1 such that the ground state of 
the system in the thermodynamic limit is an insulator 
which breaks translational symmetry. The thermody- 
namic ground state is clearly degenerate, the degeneracy 
reflecting the different broken symmetry patterns. On a 
finite lattice, such a system must thus have eigenstates 
with different crystal momentum, which, in the thermo- 
dynamic limit, become degenerate and allow us to con- 
struct linear superposition eigenstates which break the 
translational symmetry. 

Let us consider such a system on a finite lattice, with 
aspect ratio such that the thermodynamic broken sym- 
metry pattern is not frustrated. If the insulator is sta- 
bilized at a density v — p/q (with p, q having no com- 
mon factors), the flux threading argument implies the 
ground state must evolve under hc/Q flux insertion into 
a different state which has a relative crystal momentum 
AP — 2Tr{p/q)Ly, with an energy equal to the ground 
state energy in the thermodynamic limit. These states 
would be "quasi-degenerate" on a large finite lattice— 



C. Flux threading in the conventional broken 
symmetry insulator 

The manner in which the set of quasi-degenerate states 
in a broken symmetry insulator evolve under adiabatic 
flux insertion is fixed by momentum balance. Let us 
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FIG. 2: Evolution of energy levels upon flux threading in 
a conventional insulator with two-fold broken translational 
symmetry on a cylinder. The two levels which become de- 
generate ground states in the thermodynamic limit carry mo- 
menta 0, ■K. (a) The two levels cross upon threading flux along 
y in a geometry with L^ =even, Ly =odd. (b) The two levels 
return to themselves upon threading flux along y in a geom- 
etry with Lx =even, Ly =even. 



again work with a system with a two-fold broken sym- 
metry in the thermodynamic limit. 

If the filling v and Ly are such that l-nvLy = 
7r(mod 27r), flux insertion causes a momentum change 
of APa; = TT. This implies we must have two quasi- 
degenerate states differ in i-crystal momentum by tt, and 
flux threading must lead to an interchange of these two 
states. This is depicted schematically in Fig. |2Ia) where 
the two states on a finite size system, denoted by |0), Jtt), 
begin with some splitting (which must vanish in the ther- 
modynamic limit) and then evolve as the inserted flux 
$ changes. They are degenerate and cross at $ = tt 
since the Hamiltonian is invariant under time-reversal, 
but they cannot mix since the i-crystal momenta of the 
two states differ by tt even at this point. If the geometry 
is chosen such that l-cvLy = O(mod 27r) these two states 
will no longer exchange places on threading 4> = 27r (see 
Fig. [3b)). They no longer cross at $ = tt though they 
still carry relative tt momentur^^ 



D. 



Local operators to detect broken symmetry 
states 



In the presence of spontaneous translational symme- 
try breaking, there are local operators which can dis- 
tinguish the different insulating ground states obtained 
in the thermodynamic limit by taking linear combina- 
tions of the degenerate momentum eigenstates as |1) = 
(|0) + |7r))/v^ and |2) = (|0) - |7r))/V2. For example, 
broken translational symmetry along say the i, direction 
means it can be detected by some local Hermitian oper- 
ator Oi, since (l|(Oi — Oi+£)|l) ^ 0, and similarly for 
state 1 2). How does this manifest itself on a finite size 
system where such linear combinations |1),|2) are not 
eigenstates of the Hamiltonian? 



To answer this, consider the matrix element (0|Oi|7r) 
of the local operator between the eigenstates on the finite 
system. Since Oi is defined locally, it is not translation- 
ally invariant and such matrix elements will be nonzero 
in general on a finite system. However, knowing that 
T|0) = |0) and T|7r) — —\'n) we can rewrite this matrix 
element in the thermodynamic limit as 



iO.k)) ^ (l|d,|l)-(2|d,|2) 

= (i|(6,-mf-i)|i) 
= (i|(6,-o,+£)|i)^o 



(16) 



Thus, the matrix element of this local operator Oi be- 
tween states of the quasi-degenerate ground state man- 
ifold, |0) and JTr), survives in the thermodynamic limit 
and implies translational symmetry breaking. The local 
operator could be, for instance, the energy (or charge or 
current) density. 

This is the crucial difference between broken symme- 
try insulators and translationally invariant fractionalized 
insulators dealt with in the next Section. In the lat- 
ter case, matrix elements of all local operators between 
states forming the quasi-degenerate ground state man- 
ifold vanish in the thermodynamic limit. The system 
size dependence of the matrix element of local opera- 
tors between states forming the quasi-degenerate ground 
state manifold thus distinguishes an insulator with trans- 
lational symmetry breaking from a uniform fractionalized 
insulator. However, constructing such local operators 
needs some knowledge of the kind of broken symmetry, 
in contrast to our general conclusions in the earlier Sec- 
tion regarding the momenta and evolution of the quasi- 
degenerate manifold of ground states which docs not rely 
on such information. We return to this issue in Section 
EH 



In this section we have focused on conventional insu- 
lating states of a half filled system, and argued that they 
necessarily break a lattice symmetry. One case how- 
ever needs to be looked at separately, and that is the 
case of exactly tt flux through every elementary plaque- 
tte, which could be self generated in the thermodynamic 
limit. Note, this situation can preserve time reversal sym- 
metry, and hence should be admitted in our discussion. It 
is possible for such a system to be essentially an insulator 
at half filling, although it appears to possess translation 
symmetry, in that all unit cells appear identical. This is- 
sue is resolved by studying more carefully the meaning of 
translation invariance — it turns out that the operators 
that generate unit translations do not commute due to 
the presence of tt flux in the elementary plaquette. Hence, 
the smallest mutually commuting translations necessar- 
ily enclose an area equal to two unit cells, and in this 
sense we obtain unit cell doubling. Note, the emergence 
of TT flux per plaquette is a property that can be checked 
with local operators, and hence also corresponds to a 
conventional (non-fractionalized) state. 



8 



V. Z2 FRACTIONALIZED INSULATOR J* 

In the context of the insulating phase of the high tem- 
perature superconductors the question has been raised 
of whether a Mott insulator that breaks no symmetries 
could be obtained at half-filling. The analogous question 
for our bosonic system is whether a translationally invari- 
ant insulating state can be realized at half-filling. Since 
the hard-core boson state at half-filling may be viewed 
as a spin S = 1/2 system with 5*°*^^' = and C/(l) spin 
rotation invariance, this is equivalent to asking whether 
a S" = 1/2 magnet may be in a spin liquid state. The an- 
swer, after several years of work, is j/es, and one specific 
route to realizing such an insulator is via Z2 fractional- 
ization. The properties of such a phased as well as some 
microscopic models which realize them are now known^ 

The Z2 fractionalizcd insulator, X* , is a translationally 
invariant insulator. It is unconventional in that supports 
gapped fractionally charged excitations, chargons, which 
carry electromagnetic charge Qjl as well as a Z^ Ising 
charge. These chargons interact with a Z2 Ising gauge 
field in its deconfining phase. The deconfinement is re- 
flected in the presence of yet another exotic gapped neu- 
tral excitation, the Ising vortex or "vison" , which acts as 
a bundle of tt flux as seen by the chargons which carry 
Ising charge^. 

It is known that the X* phase can be realized at half 
odd-integer or integer flUings. Why doesn't the existence 
of a translationally invariant X* insulator at half-integer 
filling contradict the earlier theorem for (conventional) 
insulators? The resolution of this apparent paradox is 
that although these insulators do not break translational 
invariance, the ground state of these fractionalizcd insu- 
lators is not unique in a multiply connected geometry (in 
which the flux-threading experiment is carried out). The 
presence of the Z2 vortex, the vison, directly leads to a 
two-fold degeneracy of the ground state of the system on 
a cylinder (four- fold on a torus) . This degeneracy may be 
viewed as a result of having or not having a vison thread- 
ing each hole of the cylinder (torus). Since the vison is a 
gapped excitation in the bulk of the system, there is an 
infinite barrier for the tunneling of the vison "string" out 
of the hole of the cylinder (torus) in the thermodynamic 
limit. Thus the vison/no- vison states do not mix in the 
thermodynamic limit which is crucial to obtaining lead- 
ing a "topological degeneracy" - i.e. degeneracy which 
depends on the number of holes in the system. 

At this point we introduce the following terminology 
for the Z2 fractionalizcd insulators. The translationally 
symmetric Z2 fractionalizcd insulator at half-odd-integer 
fiUing (integer filling) will be denoted as 1*^^^ (X|j,g„). 
While in the former case, translation symmetry of a half 
filled insulator implies that the state must be exotic, it 
is of course possible to have a completely conventional 
insulator at integer filling. Nevertheless, a Z2 fraction- 
alizcd insulator may also exist at integer filling and we 
refer to this as (X*„g„). We will see below that these 
two classes of exotic insulators are in fact closely related 



to two classes of Z2 gauge theories, Z2 , ^2 



, f'^, Zl'"^" in the 
terminology of Ref . [1] 

The presence of topologically degenerate states and 
their evolution under flux threading, allows us to sat- 
isfy the momentum balance condition. The relevant case 
to consider is the translationally symmetric insulator at 
half flUing, on a cylinder with an odd number of rows. 
In this case, trivial momentum counting tells us that 27r 
flux threading leads to a degenerate state with crystal 
momentum tt. We will argue below, this momentum is 
accounted for in X*^^ since flux threading effectively adds 
a vison into the hole of the cylinder, which carries crystal 
momentum tt. 



A. Effective Hamiltonian for T* 

The effective description of X* is via a set of gapped 
charge-Q/2 bosons (chargons) also carrying an Ising 
charge, interacting with each other and minimally cou- 
pled to an Ising gauge field in its deconfining phase. In 
order to place the following discussion on a more con- 
crete footing we consider a definite Hamiltonian that can 
describe such a system, and use it to derive properties of 
the states. Since we will be interested in universal prop- 
erties that characterize the state, the results themselves 
are more general than the particular effective Hamilto- 
nian used. The simplest Hamiltonian which can describe 
a Z2 fractionalizcd insulator is: 



where 



HA{T)=Hg + H^ (17) 



H,n = -tb Y. 4 (6|6^.e-''3^-/2 + h.c.) 

m 

+ uY^{ni^2Nf (19) 



where a^'^ are Pauli matrices describing the Ising gauge 
flelds, and □ denotes the elementary plaquette on a 
square lattice. The chargons, created by b\, are min- 
imally coupled to the Ising gauge field, as well as to 
the external vector potential Aij with electromagnetic 
charge-(3/2. The second term in Hm describes repulsion 
between chargons at the same site. 

The Hamiltonian 119(1 has a local Z2 invariance under 
the transformation bi — > aiht and erf- — > aiaf^aj where 
ai = ±1. Such gauge rotations are generated by unitary 
transformations using the operator G = Y\^ Gi with 



Gi = exp [ij(l 



"»)( 



j=nn{i) 



+ 2n,)] (20) 



Local Z2 invariance implies that [Gi, iJ^(X*)] = 0. Since 
we wish to work with eigenstates of Ha{X*) which are 



invariant under such gauge transformations, translation- 
ally invariant physical states have to satisfy 



G.lphys) = (±l)|phys) 



(21) 



Let us choose Gi = +1 everywhere. 

It is instructive to first consider the limit h,U ^ K,ti,. 
In this case, since h ^ K, the gauge theory is confin- 
ing. Depending on the filling, it is then possible to show 
that one recovers conventional insulating phases such as 
a uniform band insulator (for N = even integer), or bro- 
ken symmetry states such as bond-centered (with N = 
odd integer, and U ^ h) or site-centered (with N =odd 
integer and h :^ U) charge density wave states. Thus, 
the above effective Hamiltonian in this limit is capable 
of describing well understood conventional insulators. 

However, this Hamiltonian has a richer phase diagram. 
The parameter regime where an exotic fractionalized in- 
sulator is expected for the above Hamiltonian is easily 
determined. For K ^ h, the Ising gauge field will be in 
its deconfining phase, so we can pick^^ a^^ ~ 1. Simi- 
larly, since we are interested in the insulating phase, let 
us work in the limit of large chargon repulsion U/t 3> 1 
with 2N, which is twice the filling fraction of the charge 
Q bosons, being an integer. In this limit, it is clear that it 
is energetically favorable to also set the chargon number 
Ui = 2N at each site (which is possible since 2N is an 
integer) as a starting point to understand the insulator. 
The density of charge-Q bosons in the insulator is just 
h' — N, and v could thus either be an integer or a half- 
odd integer in the X* phase corresponding to even/odd 
integer values of 27V. 

In the above regime of parameters, the system clearly 
has a charge gap 0{U) for adding an m particle (chargon) 
which is a charge-Q/2 and Ising-charged excitation that 
can propagate freely (since the gauge field is deconfined). 
It also has an energy gap 0{K) to changing a^^ — > — 1 
on a bond which changes On "'li ~* (~1) o^ adjacent 
plaquettes corresponding to creating gapped visons. It 
thus describes an exotic insulator. 



B. Flux threading in I* 

Below, we will consider the effect of threading 27r flux 
on the Z2 fractionalized insulators in the cylindrical ge- 
ometry, using the effective Hamiltonian 119|) . This will 
be done in two steps. We first consider the limit of being 
deep in the fractionalized phase (i.e. set the vison hop- 
ping to zero; h = in eauation ll9|l where it can be easily 
argued that 2tt flux threading leads to the insertion of a 
vison through the hole of the cylinder. The momentum 
balance argument then allows us to read off the crytal 
momenta of the visons in the different situations. Then, 
we turn back on a finite vison hopping /i ^ 0, and use 
continuity arguments to conclude that these crytal mo- 
menta assignments remain unchanged. 
Flux threading with static visons Consider at first 
the limit of being deep in the fractionalized phase h/K -^ 



by setting h = identically (i.e. no vison hopping), 
so that we can choose afj — 1 everywhereSS.. Let us 
adiabatically thread flux 2tt in the y direction for the 
above system on a cylinder such that the starting from 
the initial eigenstate |^(0)) in the absence of flux, the 
final state 



|*(T)) = Z^t|*(0)) 



Ut = Ttcxpi-i I HA{T,t)dt 



(22) 
(23) 



where % is the time-ordering operator. We can go to 
the Aij = gauge by making a unitary transformation 
Ha{I*,T) -^ UgHa{I* ,T)Uq'' = Ho{I*) (correspond- 
ing to zero flux). 

Since the chargons carry a charge Q/2, the hc/Q flux 
quantum threading the cylinder appears as an Aharonov- 
Bohm flux of TT for the chargons. The gauge transforma- 
tion which returns the Hamiltonian to its original form 
thus also acts on the Ising gauge flelds to remove this 
extra 7r-flux. Hence we have. 






— lAfjJJa, with 

= exp(i— ^XiUj) 






(24) 
(25) 

(26) 



and 'cut' refers to the set of links for which Xi = L^, Xj ~ 
1 (shown in Fig. 13). Thus, the flnal state in the Aij — 
gauge is |\1//) = U^UaUT\^ i) ■ Since the system is an 
insulator, the final state on threading flux <l>o rnust be one 
of the states which forms part of the degenerate ground 
state manifold in the thermodynamic limit. 

Let us deflne the loop operator Wc = Jlc '^tj '^here 
the loop C is taken around the cylinder (see Fig. |3Jl. 
Clearly, since h = Q, this operator commutes with the 
Hamiltonian (|19|l . We can use this operator to check 
whether there is a vison through the hole of the cylin- 
der. Namely, if we are in a (reference) state with erf = 1 
everywhere, then Wc = 1 and this is the no-vison state 
|u = 0). If on the other hand Wc = — 1 for each loop 
C around the cylinder, we conclude that there must be 
a vison threading the hole, and we shall refer to this as 
\v = 1). Let us evaluate Wc for the two states |^ij) 
above. We easily find W^j = 1 in the initial state. In 
order to find the eigenvalue of Wc in the final state, we 
note that since /i = 0, the initial assignment of af;: does 
not time evolve, and we only need to evaluate the effect 
of the unitary transformation Uc on Wc ■ This yields the 



/ _ 



-1. 



result that the eigenvalue in the final state is^.W^ 

Thus, for h — Q, threading a 27r flux adds a vison to 
the hole of the cylinder and interchanges the two ground 
states on the cylinder, \v = 0) ^^ \v = 1). (In fact, in the 
absence of dynamical matter fields, i.e. if, = 0, the op- 
erator Ua = riiiGcut '^ij commutes with the Hamiltonian 
and can be viewed as the "vison creation operator" intro- 
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FIG. 3: (a) Schematic figure showing a vison threading the 
hole of a cy Under in the absence of vison tunneling terms. The 
dark (light) bonds correspond to afj — —1 {afj — +1). We 
can detect the presence of the vison by evaluating the Wilson 
loop operator Yic'^lj along the contour C taken around the 
cylinder, (b) The translation operator along the x direction 
moves the dark [oi^ = —1) bonds by one lattice spacing from 
(I/x, 1) to (1,2) at each y. This is accomplished equivalently 
by acting with the gauge transformation operator d (which 
changes the sign of crfj on all bonds emanating from i) acting 
on each of the circled sites. 



ducing a vison into the hole of the cylinder and changing 
the sign of Wc-) 

Momentum balance then tells us that adding a vison 
into the hole of the cylinder must change the momentum 
of the system by 2ttv Ly. The only situation where this is 
a nontrivial crystal momentum is for the case of I*^^ on 
a cylinder with an odd number of rows (Ly odd) . Then 
we expect the two states \v = 0) and \v = 1) to differ by 
crystal momentum tt. In all other cases, i.e. for Ly even, 
or of ^*yen 7 the vison carries no momentum. 

As we shall see in the next subsection, this is consistent 
with a direct computation of the vison momentum in the 
pure Ising gauge theory. We now switch back on the vison 
hopping h ^ and ask how these conclusions might be 
affected. 

Vison state in the presence of dynamical gauge 
fields: 

Turning on a nonzero h, gives dynamics to the gauge 
field. In this case the loop product Wc no longer com- 
mutes with the Hamiltonian, we cannot use its eigenval- 
ues to label the states. Let us first see what effect this 
has on T*^^^- The two states \v = 0) and \v = 1) both 
carry zero crystal momentum, and will now mix to give 
eigenstates of the Hamiltonian. Thus, on flux threading 
there is no level crossing - threading a $o flux returns 
us to the original ground state. For X*^^ on a cylinder 
with even Ly, the two low lying states carry zero crystal 
momentum, and a similar conclusion applies. 

The situation is more interesting for T*^^ on a cylinder 



with an odd Ly. Now, the states \v — 0) and \v — 1) can- 
not mix since they carry different momenta. Thus, even 
in the presence of /i 7^ (so long as we remain in the same 
phase), we can continue to distinguish them and we can 
continue to label them as no-vison/vison states by their 
momentum, although they are not eigenstates of the Wc 
operator any longer. In this case, the crossing of the two 
states on threading a $0 flux continues to occur, since the 
crystal momentum must change by tt in order to satisfy 
momentum balance. Thus we may conclude that for the 
case of X*^^ on an odd length cylinder, the two degener- 
ate ground states (no- vison and vison through the hole of 
the cylinder) differ by crystal momentum tt. This is the 
result of the momentum balance argument applied to Z2 
fractionalized insulators. We will sometimes simply refer 
to this result for T*^^ as "the vison carrying momentum 
■K per row of the cylinder" , omitting to point out each 
time that the vison in question lives in the hole of the 
cylinder. 

The above result is consistent with the vison momen- 
tum computed using: (i) the pure Ising gauge theory (as 
shown in the next subsection), (ii) variational wavefunc- 
tions for Z2 spin liquids (as shown in Section [3D), and, 
(iii) arguments presented for short-range dimer models^ 

fori:,,. 

A side result of this analysis of identifying the vison 
crystal momenta in various situations is an unambigu- 
ous way of distinguishing fractionalized states from states 
with translation symmetry breaking for the half filled in- 
sulator. This is described in detail in Section IvH 



C. Vison momentum computed directly in the 
pure Ising gauge theory 

In order to check our deduction about the vison mo- 
mentum, let us directly compute this quantity in a pure 
Z2 Ising gauge theory without dynamical matter fields. 
If the charge gap in the insulator is large, this is the ef- 
fective description of the insulator I* . Namely, in the 
limit [/ ^ if, in the Hamiltonian H19|) and for integer val- 
ues of 2iV, a good caricature of the insulating state is to 
set Tii = N a,t each site and only consider fluctuations of 
the Ising gauge fields. This reduces the constraint on the 
physical Hilbert space to 



GT'^lphys) ^ (-l)2^exp[zj(l-a,) ^ a,^-] |phys) 



j=nn{i) 



= Iphys) 



(27) 



or equivalently, focussing only on the non-trivial case of 
cui = -1, 



n -s-=(-i) 

j=nn{i) 

in the subspace of physical states. 



2N 



(28) 
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Again, if we begin with h = 0, one ground state |u = 0) 
of the gauge theory on a cyhnder may be obtained as 
the reference state a^ — 1 projected into the physical 
subspace, and for this one has the loop operator Wc — 1 • 
A second (degenerate) state may be obtained by acting 
on this ground state with 



V , = 



n 






(29) 



which commutes with the Hg for h = 0. The subscripts 
{Lx, 1) on yt are a mnemonic for the column on which 
the a^ operators act as shown in Fig. |2Ia) . The resulting 
state has Wq = — 1- Let us compute the momentum of 
these two states. Clearly, the state \v = 0) has zero mo- 
mentum since it is translationally invariant by construc- 
tion. To compute the momentum of the second state, we 
first note (see Fig.O^b)) that 



TVL.iT 



V12 



m n 



i.i-\-x ) 



jenn{i) 



V, 



ix,l 



(30) 



(31) 



where the sites i have Xi = L^ and correspond to the 
circled sites in Fig. ^h) . Using the constraint in Eq. |2S1 
this reduces to 






n( 1: 



2N 



y/^ 1 = exp{i2TTNLy)Vl^j^. 

(32) 
Thus, for the second state, \v ^ 1) = V^ -^|w = 0), acting 
with the translation operator leads to 



T\v^l) = 



4 



TVL.iT- 



= exp{i2TTNLy)vl2\v = 0) 



ex.p{i2TTNLy)\v = 1) 



(33) 

(34) 
(35) 



In other words, (i) for even 2N, namely in I*,je„, the 
state |f = 1) carries zero crystal momentum and , (ii) for 
odd 2iV, namely in X*^^, the vison state \v = 1) carries 
momentum irLy. 

As before, we can now turn on a nonzero field h. In 
I*j,g„ , the two ground states will mix and split in a finite 
system, since they carry the same momentum quantum 
number. The same is true for T*^^ with even Ly. How- 
ever, for X*^^ with odd Ly, the two ground states carry 
relative momentum tt, thus they cannot mix even on a 
finite system with nonzero h and can be distinguished by 
their momentum. 

Finally, we may introduce dynamical matter fields. Al- 
though the operator in Eqn. 1291 no longer can be iden- 
tified as a vison creation operator, the low energy struc- 
ture of the system, i.e. topological degeneracies will not 
change as long as we are in the same phase. Also, since 
the crystal momenta of these low lying states on the 
cylinder can only be one of 0,7r (from time reversal in- 
variance), continuity requires that the crystal momentum 



assignments made before for the low lying states continue 
to hold in the presence of dynamical matter fields. This 
constitutes a direct check of the results deduced above 
using momentum balance arguments. 



D. Momentum computation from variational 
wavefunctions for the vison 



So far, we have discussed bosonic models for the insu- 
lators X* . However, as mentioned earlier, we can view 
a hard-core boson as a S* = 1/2 spin, and the insulat- 
ing state I*^^ with 27V an odd integer as a Z2 fraction- 
alized spin liquid insulator. Such spin liquid insulators 
have long been of interest in connection with frustrated 
magnets and the high temperature superconductors. The 
Q — 1/2 chargon excitations in the bosonic language cor- 
respond to 5 = 1/2 excitations (called spinons) in the 
spin liquid. What do the visons in 1*^^ correspond to? 

To answer this, we note, following AndersonS^, that 
one can represent of the ground state wavefunction for 
spin liquids with short-range antiferromagnetic correla- 
tions by "Gutzwiller projecting" a superconducting wave- 
function, i.e. restricting to configurations with a fixed 
number of electrons per site. Such a picture also emerges 
from mean-field studies of frustrated magnets using a 
fermionic representation for the spins. This suggests that 
perhaps excitations of the spin liquid may also be related 
to excitations in the superconductor. Following this line 
of thought, the S = 1/2 spinon in the spin liquid may 
be viewed as a projected Bogoliubov quasiparticle of the 
superconductor. Similarly it is natural to expect that the 
hc/2e vortex in the superconductor becomes the vison2&. 

We can check this possibility by computing the mo- 
mentum of a projected hc/2e BCS vortex threading the 
cylinder, with odd/even number of electrons at each site, 
and seeing if it agrees with the results for l^^.d/'^even ob- 
tained above. To do this, we write the BCS state (with 
total electron number Nf.) as 



\BCS{N,))=[Y,. 



iVe/2 



'k^Tcika 



(36) 



where 0k denotes the internal pair state of the Cooper 
pair formed by (k, "f ) and (— k, |). which carries zero 
center of mass momentum. 

To get the spin liquid state at half-filling, we have 
to choose Ne = N, the number of lattice sites, and 
Gutzwiller project this state by acting with the opera- 
tor Pg — Yiii^ ~ '^iT'^u) which eliminates configurations 
in which two electrons occupy the same site. The vari- 
ational ansatz for the spin liquid ground state is then 
Pg\BCS{N)), and it is a translationally invariant state 
with zero momentum. To construct a BCS hc/2e vor- 
tex threading the cylinder, we need to pair states with 
(k + q/2, t) and (-k + q/2, i) with q = {2n/L^)x and 
amplitude (f>u_, this leads to the iV-particle vortex state 
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\hcl2e{N)) carrying a momentum of q per pair, or a to- 
tal momentum (21:/ Lx){N /2). Setting N — L^Ly, we 
see that the momentum of this state is just 'KLy. The 
projection operator Pq commutes with the translation 
operator. Thus the trial vison state \v) — PG\hc/2e{N)) 
has a momentum irLy in agreement with earlier argu- 
ments for the X*^^ insulator. 

With an even number of electrons at each site, the 
vison wavefunction carries no crystal momentum. This 
is consistent with our earlier result for J* ,„„ . 



VI. IDENTIFYING Z2 FRACTIONALIZED 
STATES IN NUMERICAL EXPERIMENTS. 



Numerical investigation of microscopic models, for ex- 
ample exact diagonalization studies^, are an important 
tool in finding new states of matter such as states with Z2 
topological order. In this context it is important that reli- 
able diagnostics be available for the identification of these 
states in the system sizes that can currently be solved 
on the computer. This question may seem straightfor- 
ward in principle, the four fold topological degeneracy 
of the Z2 states on a torus that are indistinguishable by 
any local operator seem to provide a unique prescription. 
However, in practice there are several potential problems. 
First, since the numerical simulations are performed on 
finite sized systems, states that are degenerate in the 
thermodynamic limit are only approximately so in these 
systems. The problem is particularly severe when gap- 
less Z2 charged matter fields are present, in which case 
the splitting between topological sectors that differ by 
the presence of a vison can be large and go down to zero 
only algebraically with system size (in contrast to the ex- 
ponentially small splitting in the absence of such gapless 
gauge charged matter fields). Secondly, if degeneracies 
arise as a result of broken translation symmetry, rather 
than topological order, the relevant order parameter for 
this translation symmetry breaking may be hard to iden- 
tify, and hence we would like to have available a pre- 
scription for distinguishing such states even if the order 
parameter is not known. Below we use insights from the 
momentum balance arguments to resolve both of these 
issues. Indeed we will see that the analysis of the previ- 
ous sections, with their focus on finite sized systems and 
crystal momentum quantum numbers, are ideally suited 
to addressing these questions. We begin by addressing 
the second of these two questions first - i.e. given a set of 
states comprising the low energy manifold of the system 
as the thermodynamic limit is approached, how does one 
distinguish topologically ordered states from a conven- 
tional translation symmetry broken state? 



A. Z2 Topological Order vs. Translation Symmetry 
Breaking 

Consider a system of bosons on a lattice at half filling 
(or equivalently a spin 1/2 system with one spin per unit 
cell) . As discussed previously, a translationally invariant 
insulating phase implies the presence of topological or- 
der (although it is possible to have topologically ordered 
phases that also break translation symmetry). Assume 
that a group of low lying states have been identified — 
under what conditions can we associate these with the 
degeneracies assocated the Z2 topological order, rather 
than with low lying states leading to translation symme- 
try breaking? 

First consider the system in the cylindrical geometry 
with an odd number of rows {Ly odd, L,j: even in Fig.^l. 
Then, flux threading ensures that we will have two low ly- 
ing states with crystal momentum P^ ^ 0, Px ~ it which 
are interchanged on threading 27r flux. This is irrespec- 
tive of whether the system is heading towards translation 
symmetry breaking or towards a Z2 topologically ordered 
state in the thermodynamic limit. Therefore this setup 
is not particularly useful for discussing for distinguish- 
ing the two states. One may also consider the toroidal 
geometry for an even x odd system, however, this can 
potentially frustrate certain patterns of translation sym- 
metry breaking in a half filled system and hence we do 
not consider it further here. 

Now consider the system in a toroidal geometry, but 
with both Ly, L^ even. Now, we have shown earlier that 
a Z2 topologically ordered state will have four low ly- 
ing excitations, all with zero crystal momentum in this 
geometry. A conventional translation symmetry break- 
ing state on the other hand will invariably have at least 
one state in the low energy manifold that carries nonzero 
crystal momentum, in addition to a zero crystal momen- 
tum state. This follows almost by definition, in order 
to build up a translation symmetry breaking state one 
needs to make a linear combination of states with differ- 
ent crystal momenta. This then is a precise way to tell 
apart a Z2 topological state from a more conventional 
translation symmetry breaking state, which just requires 
a correct identification of the low energy manifold and 
the crystal momenta of these states. The translation- 
ally symmetric topologically ordered state is present if 
there are four low lying states with zero crystal momenta. 
If topological order coexists with translation symmetry 
breaking, then too this quadruplet of zero momentum 
states persist, although other quadrupled states with dif- 
ferent crystal momenta will be present in the low energy 
manifold. This is true for both I^ven ^^^ ^odd- 
Recent exact diagonalization studies studies of a mul- 
tiple spin exchange model on a triangular lattice have 
found signatures of an interesting new spin state, which 
has been proposed to be a topologically ordered spin liq- 
uid phase in a certain regime of parameters^. Let us ap- 
ply the method of distinguishing topological order from 
broken translational symmetry discussed above, to these 
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states. 

In the parameter regime of interest, the system in 
Ref. was argued to be heading, in the thermodynamic 
limit, towards a spin gapped phase without long-range 
magnetic order. Furthermore, a set of three spin sin- 
glet states which appear to become degenerate with the 
ground state with increasing system size were identified. 
The authors were unable to find simple valence bond 
(eg. nearest neighbour) crystal states that would lead 
to degenerate ground states with the quantum numbers 
(crystal momenta, rotations, reflections) of these low ly- 
ing states. Hence they identified this apparent four-fold 
degeneracy with the degeneracy arising from topological 
order of a Z^ fractionalized spin liquid, such as described 
by Eqn^l on a torus. While this would be a very in- 
teresting result, we can ask if the quantum numbers of 
these nearly degenerate states are consistent with those 
of a vison in an odd Ising gauge theory, that we have 
derived earlier. However, the three excited states which 
appear to become degenerate with the ground state carry 
nonzero crystal momentum on an 6 x 6 lattice^. This 
is in disagreement with our conclusion regarding vison 
states in a Zi fractionalized phassSI^ namely, that they 
carry zero crystal momentum on evenxeven lattices. We 
therefore conclude that this interesting identification of 
Z2 topological order in these systems does not stand up 
to detailed scrutiny. The actual nature of the phase being 
approached by these systems then remains an open ques- 
tion, especially since an extensive search of conventional 
broken symmetry states in Ref _4 did not yield a candidate 
phase. The remaining possibilities are perhaps a con- 
ventional translational symmetry broken valence bond 
crystal phase, involving non-nearest neighbour dimcrs, 
some other more exotic fractionalized state, or that the 
all the low lying states associated with the broken sym- 
metry have not been identified as a consequence of finite 
size effects. Note, the evolution of these states under 
flux threading which they have studied on odd x even and 
evenxeven lattices is also consistent with a conventional 
broken symmetry state. 



B. Eliminating the Spinon Contribution to Vison 
Splitting 

In this subsection we will utilize the flux threading pro- 
ceedure to find a way of eliminating the splitting between 
the topologically 'degenerate' states that arises from the 
presence of Z2 charged matter fields. While in the ther- 
modynamic limit, a Z2 fractionalized system must posess 
a four fold ground state degeneracy on the torus and a 
two fold ground state degeneracy on the cylinder, in a 
finite system this splitting may be so large that it makes 
the identification of the low energy manifold problem- 
atic. In this subsection we will utilize the fiux threading 
proceedure to find a way of eliminating the part of the 
splitting between these states that arises from the pres- 
ence of Z2 charged matter fields. 



In order to study the problem in more detail con- 
sider a finite sized system in the cylindrical geometry 
that is heading towards a Z2 fractionalized insulating 
state. Consider first the situation on an evenxeven 
lattice. The low energy manifold consists of a pair of 
states that eventually become degenerate in the thermo- 
dynamic limit, but at this stage have a finite splitting 
Ai?. The spitting arises from two sources: first there is 
vison tunneling, that mixes the zero and one vison states, 
which acting alone, would lead to a splitting of Ai^hop- 
Since this involves a gapped vison (with gap e) hopping 
across the entire height Ly of the cylinder, one would 
expect this to be exponentially small in their product 
i.e. Aii^hop oc exp(— cLj,e), where c is a constant. The 
second contribution to the splitting arises from the pres- 
ence of matter fields that carry gauge charge. Clearly, 
the presence or absence of a vison will affect the propa- 
gation of these particles and in the absence of vison tun- 
neling will give rise to an energy splitting Ai^mat- With 
gapped Z2 gauge charged matter fields (with a gap e), 
clearly this splitting will require virtual processes where 
the gapped particle goes once around the cylinder which 
implies Ai^mat oc exp(— c'L2;e). The total splitting is 

easily seen to be: AS = J l\E\^^ -K A£;^af Thus, in 
situations like the one described above, where both vi- 
sons and gauge charged matter have a healthy gap, finite 
sized system studies can in practice isolate the low en- 
ergy multiplet that leads to topological degeneracy in the 
thermodynamic limit iS Note that the presence of gapless 
matter fields which are gauge neutral do not affect these 
conclusions. Further, atleast in principle, the topologi- 
cal degeneracies described here, with splittings which are 
exponentially small in system size can be separated from 
low-lying modes of the gauge neutral excitations whose 
splitting scales inversely with system size. 

However, if the gauge charged matter is gapless (e.g. if 
there are fermionic Z2 gauge charged excitations with a 
Dirac spectrum that often appear in mean field theories 
of spin liquids) then the splitting of the zero and one vi- 
son states are no longer exponential in the perimeter size 
of the system, but only a power law i.e. A_Einat oc L~'^ , 
where 77 > 0, and this dominates A£'. This is poten- 
tially a serious problem since in a finite sized system the 
splitting is very likely to be large and and also hard to 
distinguish from the low energy states arising from the 
gapless fermions, which also have energies that vanish 
as the inverse size of the system. Below we will prove 
that in the presence of vr flux (i.e. antiperiodic boundary 
conditions for the unfractionalized bosons/magnons) the 
matter contribution to the vison splitting is switched off! 
Essentially this arises because the Z2 charged particles, 
which also carry half a unit of charge, see the antiperi- 
odic boundary conditions as a flux of (say) — 7r/2. Adding 
a vison then implies a flux of 7r/2. However, since these 
two situations are related by time reversal symmetry, the 
energy contribution from the matter fields in these two 
cases is identical - which implies that the splitting arises 
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solely due to vison tunneling, which can be made expo- 
nentially small. 

In order to show that at a flux of tt the matter contribu- 
tion to the vison splitting vanishes Ai^mat = 0, we adopt 
the following procedure. We consider the system (head- 
ing towards a Z2 fractionalized state) on an evenxeven 
cylinder, and consider first the limit where the vison hop- 
ping is turned off. Then, the splitting of the low energy 
states occurs entirely because of the gauge charged mat- 
ter contribution. We now consider introducing a vison 
through the hole of the cylinder and argue below that ex- 
actly at flux TT the two states are exactly degenerate (even 
in a finite system). Since the vison hopping has been 
tuned to zero, the remaining source of splitting (arising 
from the gauge charged matter fields) must also be zero 
at this point. We can then reintroduce the vison hop- 
ping, and the states at flux tt will now be split, but the 
splitting occurs entirely from vison tunneling. 

Let us now show that in the absence of vison tunnel- 
ing, an exact degeneracy occurs at tt flux. We consider 
for deflniteness the model in Ean. H19|) - although it con- 
tains gapped matter flelds with Z2 gauge charge, the con- 
clusions simply show that all matter field contributions 
cancel at this special fiux, and hence can be easily ex- 
tended to the case of gapless matter fields as well. We 
consider the limit of vanishing vison tunneling (i.e. h — Q 
in Eqn. I19|) . We start with zero flux through the cylinder 
and consider inserting, adiabatically, a flux of I-k. 

Then, as argued in Subsection inserting a 27r flux 
leads to insertion of a vison. Now, at zero flux, the two 
low lying states can be classified in terms of vison num- 
ber, since the vison hopping has been set to zero. The 
vison number is measured by the operator Wc = Jlc '-'I] 
where the loop C is taken around the cylinder (see Fig. 2). 
Clearly, since /i = 0, this operator commutes with the 
Hamiltonian, and the two low lying states can be labelled 
with the eigenvalues of (1 — Wc)/2, i.e., the vison num- 
ber. The splitting between these levels arises entirely 
from the gauge charged matter fields. On flux threading, 
these two states must then interchange - since threading 
27r flux inserts a vison in this limit. This means that the 
two levels have to cross at some point (or more gener- 
ally at an odd number of points) as a function of flux. 
Now, time reversal symmetry tells us that if there is a 
crossing point at flux 0, then there must be one also at 
the point 27r — 0. Thus, in order to arrange for an odd 
number of crossing to ensure the levels do interchange, 
we need that there is always a crossing at flux tt. Thus, 
the two states with vison and no vison are exactly de- 
generate at this value of the flux and hence we conclude 
that the splitting from the gauge charged matter flelds 
vanishes at this value of flux. Now, turning on the vison 
hopping h ^ Q will lead to a flnitc splitting even at tt flux, 
but this splitting arises entirely from vison tunneling and 
hence at this value of the flux we have Ai? — Ai^hop and 
hence vanishes exponentially in the width of the system. 
These arguments can easily be taken over to the toroidal 
geometry as well. 



We note that this result is useful even in the study 
of SU(2) symmetric spin liquid statesi, where although 
the intoducing the tt flux will require breaking the SU(2) 
symmetry, this only occurs along one row of the cylinder 
(e.g. changing the sign of the exchange constants for the 
SxSx and SySy interactions), and hence may be viewed 
as a fairly weak perturbation away from full SU(2). It 
should also be useful in projected wavefuntion studies, es- 
pecially in establishing the existence of Z2 fractionalized 
states with excitations that have a Dirac dispersion2&. 

Finally we note that while turning on a flux of tt is ef- 
fective in cancelling the splitting arising from dynamical 
matter flelds, the splitting from vison tunneling can be 
cancelled in a like manner by considering a cylinder with 
an odd number of rows at half fllling, where the vison 
and no vison states differ by crystal momentum tt and 
hence do not mix. When both these processes are active 
we expect the vison and no-vison states to be exactly 
degenerate. Indeed, this is borne out by the observation 
that for a cylinder with an odd number of rows at half 
fllling, when the threaded flux reaches tt there is always 
a level crossing just from momentum balance arguments 
and time reversal symmetry (see Fig. 2(a)). 



VII. MOMENTUM BALANCE FOR FERMI 
LIQUIDS 

A. Conventional Fermi Liquids 

Let us first briefly review the momentum balance ar- 
gument due to Oshikawa^ for conventional Fermi liq- 
uids where it leads to Luttinger's theorem. Consider 
fermions with charge Q and spin t, | at a fllling per site 
of v^ 



n 



u. Now consider flux threading in the 



cylindrical geometry of Fig. Q] with L^ columns and Ly 
rows. We imagine threading unit flux $0 — ft-c/Qthat 
only couples to the ] spin fermions. Via trivial momen- 
tum counting this proceedure can be seen to impart a 
crystal momentum of: 



APj; = 2'KV^Ly 



(37) 



Similarly, one could imagine performing the flux thread- 
ing with the cylinder wrapped along the perpendicular 
direction which would yield a crystal momentum change: 



APj, ~ 2tiv^Lx 



(38) 



Now in the regular Fermi liquid phase, this crystal mo- 
mentum imparted during flux threading is accounted for 
entirely by quasiparticle excitations that are generated 
near the Fermi surface. Using the fact that long lived 
quasiparticles exist near the Fermi surface, and the fact 
that the Fermi liquid is adiabatically connected to the 
free Fermi gas, the quasiparticle population 5np excited 
during the flux threading proceedure can be worked out. 
Clearly, flux threading for noninteracting fermions will 
lead to a uniform shift of the Fermi sea by Ap^; = 2tt/Lx 
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from which the quasiparticle distribution function can 
be determined. Indeed all of these excitations are close 
to the Fermi surface, which is required in order to apply 
Fermi liquid theory. The total crystal momentum carried 
by these excitations can be written as: 



A> = ^S 



UpP 



(39) 



It is convenient to first evaluate this expression neglecting 
the dicrete nature of allowed momentum states in a finite 
volume system and treating the shift in the Fermi sea 
dpx — 27r/La; as infinitesimal. This yields: 



from unity). Then, m^Lx = rUyLy 
multiples of L^Ly] namely rUxL^ = 



implies that they are 
pL 



lyLy 



p an integer. Thus we obtain the result: 



Ly with 



Vfs 

{2^f 



P 



(44) 



which of course is Luttinger's theorem- that relates the 
Fermi volume to the filling (modulo filled bands that are 
represented by the integer p) . 



VIII. MOMENTUM BALANCE IN FL* 



APr 



Sp ■ dSp 



(40) 



where dSp is a vector normal to the Fermi Surface, and 
the integral is taken around the Fermi surface. Using 
Gauss divergence theorem this can be converted into an 
integral over the Fermi volume which yields: 



AR 



Spa 



dV 

27r 27r 

FV T-TZ 



(41) 



thus, 



APx = ^^ 



j^ 271- 2-n- 



A more careful derivation that keeps track of the dis- 
creteness of the allowed momenta gives the same result. 
Clearly the relevant Fermi volume that enters here is that 
of the up spins. Below we assume for simplicity that both 
the t and [ spins are at equal filling and so V] — vi — v 
and V^g = Vpg = Vps- Equating the results from the 
trivial momentum counting, and the momentum count- 
ing above for the Fermi liquid (up to a reciprocal lattice 
vector) yields: 



2'Ki'L 
2ttvL 



27r Vfs 



2^ Vfs 



Ly (27r) 



LxL,, + 2-Kmx 
^LxLy -I- 2'iTmy 



where m^ and ruy are integers and the two equations 
above are obtained from threading flux in the x and y 
directions. These equations can be rewritten as: 



N - LxLy = LxTTlx 

N - LxLy—-— = Lyniy 



(42) 
(43) 



where we have introduced the particle number N — 
LxLyV, an integer. In order to obtain the strongest con- 
straint from these equations we consider a system with 
Lx , Ly mutually prime integers (no common factor apart 



Here we will consider an exotic variant of the Fermi 
liquid, where electron like quasiparticles coexist with Z2 
fractionalizationSi. This state may be obtained begin- 
ning with a Z2 fractionalized insulating state of electrons 
that breaks no lattice symmetries. We consider a specific 
model where the spinous (/J spin half, charge neutral ex- 
citations) are fermionic and the chargons (6^ spin zero, 
unit charged excitations) are bosonic. The electron oper- 
ator is written as cj. = bf^ and the relevant gauge struc- 
ture is Z2 which implies that pairing of spinous is present. 
If one is deep in the insulating phase then there is a large 
gap to the chargons; furthermore, if there is also a spin 
gap, then the low energy effective theory is just an Ising 
gauge theory. For an insulator with an odd number of 
electrons per site, in this regime we may set the chargon 
number Uf, — and the spinon number nf = I. The Ising 
gauge charge at each site is (—1)"''+"/, which leads to an 
odd Ising gauge theory in this situation. For an insula- 
tor with even number of electrons per site, an even Ising 
gauge theory would result. We now imagine a situation 
where the lowest charge carrying excitation in the system 
is the electron itself. This could arise if the spinon and 
chargon form a tightly bound state so that it has a lower 
net energy than an isolated chargon. Doping would then 
lead to a 'Fermi liquid' of electron like quasiparticles, co- 
existing with gapped visons, spinous and chargons, which 
is the FL* phase we wish to discuss. It already appears 
that a violation of Luttinger's relation may be expected 
here if we dope an insulator with an odd number of elec- 
trons per unit cell, since only the doped electrons may 
be expected to enter the Fermi volume. Here we will 
see how momentum balance arguments allows for such a 
violation, but nevertheless constrains the possible Fermi 
surface volumes so that a generalization of Luttinger's 
theorem to this exotic class of Fermi liquids holds. 

In order to follow in detail the evolution of the sys- 
tem under flux threading we study the following model 
Hamiltonian: 



Hfl' -- 


- iJg + H^p_^f^ + Hint + Hgauge (45) 


Hi ~ 


' E ^?.4^^- (46) 




<ij> 




= -E*?;4&l^.-E*?;4/t/.- 




<ij> <ij> 
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H, 



gauge 



i 



with the constramt on all physical states: 

j=nn{i) 



(47) 
(48) 

(49) 



U^ 



and Hint denotes the interactions between various which 
we do not specify here except for assuming that terms 
here do not couple to an externally imposed gauge field 
that is required for flux threading. Note, the spinous 
and chargons are coupled to the Z2 gauge field, while 
the electrons are of course Z2 gauge neutral. We have 
selected, for simplicity, an on site pairing interaction for 
the spinous; while such onsite pairing terms are absent 
in microscopic models that forbid double occupancy of 
electrons, here we are concerned with universal aspects 
of quantum phases which are not affected by this simpli- 
fication. 

Flux Threading in FL*: We now consider the effect 
of flux threading on the ground state of the Hamilto- 
nian (|48|l . In the presence of a vector potential A^ {A^ ) 
coupling to the up (down) spin electrons, the hopping 
matrix elements for the up (down) spin electrons is mod- 
ified to tfj -^ tfjC^ *3 {tf^ — * tfjC^ 'j) , while the chargon 

hopping amlitude is modified to i^ — + t'^.e'^^ ^^ ^i' and 
the up (down) spinon hopping amplitude is modified to 

t% -^ ifj-e5(<-^»'.) [tlj -> tf^.e-5(<-^>'.)). Below we 
imagine threading 2tt flux in A^ and study the evolution 
of the ground state of the system in this process. In ad- 
dition to the excitation of particle-hole pairs of the elec- 
tron like Fermi-liquid quasiparticles, we will also see that 
in some situations a vison excitation is inserted through 
the cylinder which gives rise to the modifled Luttinger 
relations. We begin by considering flux threading in the 
absence of vison dynamics (ft, = in Eqn. I48|) . where 
results are easily derived, and then reinstate the vison 
dnamics and show that the central result is unaffected. 

The adiabatic insertion of a unit flux quantum that 
couples to the up spin electrons is affected by introducing 
a gauge field on the horizontal links of the cylinder in 
FigureHand increasing its strength from zero (A| = ^ 
27r/ij;) in time T. The time evolution of the quantum 
state can be written as \'4){T)) = Z/^r|'0(O)) where Ut = 

Tt exp ( —i /g HpL* {t)dt j where Tt is the time-ordering 

operator, and the time dependence of the Hamiltonian 
arises from the flux threading. Clearly, since 27r of flux 
is invisible to the electrons, the final state must be some 
excited state of the initial {A^ = 0) Hamiltonian. In 
order to make this explicit, the 27r fiux is gauged away, 
which can be accomplished by the operators UaU^ with 



ij ^cut 



(50) 



U^ = exp 



*7^E^*K«+ 9(' 



'/T 



Vi 



<)) 



While the first unitary operator eliminates the gauge 
field for the electrons, it changes the sign of the hopping 
matrix element on a single column of horizontal links for 
the chargons and spinous which behave like half charges. 
This modification to the hopping can be absorbed in the 
Z2 gauge fields, which is accomplished by the unitary 
operator U^, which returns us to the initial Hamiltonian. 

The action of the time evolution operator Ut is to ex- 
cite electron like quasiparticles about the Fermi surface 
in the usual manner, while the gapped spinous and char- 
gons are not excited during this adiabatic flux threading. 
In the absence of vison dynamics {h — in equation I48|l. 
it may be seen that a vison is also introduced during the 
flux threading proceedure. This is argued as follows. In 
the absence of vison dynamics, the vison number through 
the hole of the cylinder, as measured by the operator 
Wc — Tlc'^ij, where C is a contour that winds around 
the cylinder, is a good quantum number since it com- 
mutes with the Hamiltonian in this limit. However, in 
the course of flux threading and returning to the original 
gauge, it changes sign since it may be easily veriflcd that 
UaWcU~^ — —Wc, which implies vison insertion. Thus, 
the final state has a displaced Fermi sea and a vison. 

We can now combine the results of trivial momentum 
counting and a knowledge of the vison momentum to ob- 
tain the volume of electron like quasiparticles. This is 
most easily done in the limit of a very large gap to the 
gauge charged particles (spinous and chargons). Then, 
the phase is described by gauge neutral electrons form- 
ing a Fermi liquid and a pure Ising gauge theory in the 
deconfined phase. We know that the vison excitations 
of the latter through the hole of the cylinder with an 
odd number of rows, carries crystal momentum or tt 
depending on the even or odd nature of the Ising gauge 
theory. The gauge constraint in Eqn. 0^ tells us this 
depends on the parity of Uf -f rii, at each site. If we 
set n\, — Q for the gapped chargons and Uf — 1 for the 
gapped spinous, where for the latter we assume the sys- 
tem is obtained continuously by doping a spin liquid with 
spin 1/2 per unit cell (i.e. a spin version of 7*^^^). In this 
limit an odd Ising gauge theory will be obtained , where 
the vison threading an odd width cylinder carries crystal 
momentum tt. If the gap to the spinous and chargons 
is now reduced from infinity, this crystal momentum as- 
signment to the vison cannot change continuously (from 
time reversal symmetry) , and is hence expected to be in- 
variant for a finite range of gap values. Thus, the phase is 
expected to be continuously connected to the large gap 
situation with integer or half integer filling, which will 
determine the momentum assignments. 

Thus, we are left with the result that two types of ex- 
otic Fermi liquid states f L*„g„ and FL*^^ are expected, 
that differ in the crystal momentum carried by the vi- 
son excitations. The momentum balance argument then 
immediately implies that these two states will have dif- 
ferent Fermi volumes at the same filling — while FL^^^^ 
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will have a Fermi volume that is identical to that of a 
conventional Fermi liquid at the same filling and hence 
respects Luttinger's relation, FL*^^^ has a Fermi volume 
that violates Luttinger's relation in a very definite way. 
Since the only situation where the momentum balance 
argument will give a result distinct from that of a con- 
ventional Fermi liquid is for the case of FL*^^^ on an 
even X odd lattice, where the vison carries a nontrivial 
crystal momentum, we discuss that below. Consider fiux 
threading in such a phase on a cylinder with an odd num- 
ber of rows. This will introduce a vison through the hole 
of a cylinder which carries crystal momentum tt. This 
needs to be subtracted from the usual momentum bal- 
ance relations for a Fermi liquid displayed in Eqn. 1431 
This leads to a modified Luttinger relation between the 
Fermi volume in FL*^^^ and and the electron filling v. 



1, 
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(51) 



where p is an integer that represents filled bands. The 
crucial difference from the usual Luttinger relation in 
Eqn. ^J is the fact that the Fermi volume is determined 
by i^ — i , which is related to the fact that it is obtained 
by doping the fractionalized spin model which is transla- 
tionally symmetric at half filling. 

Finally we argue that reintroducing the vison hopping 
[h 7^ 0) does not affect these conclusions. In the cases 
where the vison threading the cylinder carries zero crystal 
momentum, introducing vison hopping leads to a mixing 
of the vison and no vison states in a finite system. This 
implies that we are no longer guaranteed to have a vi- 
son on flux threading. However, for the case of FL*^^^ 
on a cylinder with an odd number of rows, where the 
vison carries crystal momentum tt, the nontrivial crystal 
momentum blocks the tunneling of visons even on a fi- 
nite sized system. This implies that flux threading does 
indeed introduce a vison which finally leads to the mod- 
ified Luttinger relation. We can see that the tt crystal 
momentum carried by the vison cannot be transferred 
to the only other gapless excitations in the problem, the 
electron like quasiparticles, since they are gauge neutral. 



IX. CONVENTIONAL SUPERFLUID ST 

Luttinger's theorem was formulated for Fermi liquids, 
and we have extended Oshikawa's argument to show how 
the theorem must be modified to account for the ex- 
istence of gapped spin liquid insulators (which may be 
equivalently viewed as fractionalized bosonic insulators) 
as well as fractionalized Fermi liquids. In both cases, the 
presence of topological order was crucial. Let us next 
turn to conventional superfluids and ask: What property 
of the superfluid phase is captured by the Oshikawa ar- 
gument, and gets fixed by the particle density? While we 
focus on the case of bosonic superfluids, we expect our re- 
sults to also be applicable to s-wave superconductors with 



a large gap, so that the resulting Cooper pairs may be ef- 
fectively viewed as bosons. Also, we consider the case of 
neutral bosons (no internal electromagnetic gauge field). 
Again, these results could be applied approximately to 
the case of charged superfluids if the penetration depth 
is sufficiently large. 

Conventional superfluids in two or more dimensions 
[D > 2) are Bose condensed at zero temperature, and 
have a unique ground state (on both cylinders and torii). 
It is clear that the Oshikawa argument must then cap- 
ture some property of the excitations in the superfluid. 
A conventional superfluid supports two kinds of excita- 
tions: the gapless linearly dispersing Goldstone mode of 
the broken symmetry state ("phonon"), and topological 
defects, namely vortices. We show below that it is the 
Berry phase acquired by a vortex on adiabatically going 
around a closed loop that is fixed by the particle density 
v, independent of the strength and nature of interactions 
between bosons in the superfluid. 



A. Effective Hamiltonian for SJ- 

We may describe a conventional superfluid most con- 
veniently in a rotor representation for the bosons 



thus Bl 



-i(f)j: 



, BlBr -> Nr with [exp{i(t)r),Nr'] 



exp{i(j)r)Srr' , and the Hamiltonian for interacting bosons 
in these variables takes the form 






cl>j+QA,j) + Vint[n] (52) 



where the interactions may be of the general form 
14nt[?^] — X^rr' ^rr'^r^r'- The Bose coudenscd super- 
fluid, which results in dimensions D > 2 when tb is the 
largest scale in the Hamiltonian, supports linearly dis- 
persing "phonons" , which are the Goldstone mode of the 
broken symmetry. Vortices appear as topological defects 
in the phase field in the ordered state, and there is a 
nonzero gap to creating vortices in the bulk of the super- 
fiuid. 

In the above discussion, we have assumed that the 
phase variable has periodic boundary conditions, namely, 
ifir+L^ = 'Pr, 'Pr+Ly = fr- Howevcr, On cylindcrs/torii 
the superfluid has additional excited states corresponding 
to creating vortices through holes of the cylinders/torii. 
A state with (fr+L^ — fr + 2'Kmx corresponds to a 
strength vrix vortex through a hole in the cylinder. 



B. Flux threading in ST 

Consider N bosons each with charge Q condensed into 
a conventional superfluid ground state on an L^ x Ly lat- 
tice in the form of a torus. Trivial momentum counting 
tells us that threading flux $o = hc/Q into the cylinder 
on which the system lives changes the crystal momentum 
by 2-KvLy^ where u is the filling and Ly is the number 
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of rows of the cylinder. We show below, using a low- 
energy description of the superfluid, that adiabatic flux 
threading introduces a vortex into the hole of the cylin- 
der/torus. We do this two steps. First, we turn off the 
boson interactions which allows us to directly construct 
the final state and see that it corresponds to introducing 
one vortex. Next, we turn back on the boson interactions 
and argue that this does not affect the state or change 
the momentum carried by the vortex. 

Threading flux $o introduces a vortex: 

Let us adiabatically thread flux hc/Q in the —y di- 
rection for the above system on a cylinder as in Fig. ^ 
such that starting from the initial eigenstate |^E'(0)) in the 
absence of flux, the final state reached is \^{T)). The fi- 
nal state can of course be written as \^{T)) = Ut\^{Q)) 

with Ut — 7^ exp ( — z Jg HA{SJ-,t)dt] where % is the 
time-ordering operator. We can go to the Aij = 
gauge by making a unitary transformation Hj^ {SJ-, T) — > 
UgHa{SJ^,T)Uq^ = Hq{ST) (corresponding to zero 
vector potential). Here IAg = l^4>, with 



Kd, 



exp(«— y ^XjUi 



(53) 



Thus, the final state in the Aij = gauge is |^/) = 

Let us consider the extreme limit achieved by l^nt — *■ 
0. In this case, the system is initially in the full Bose 
condensed state \exp(i(f)i) — 1), which is unaffected by 
the time development operator Ut, and the final state 
after acting with Utj, has \exp(i(f)i) — exTp{i2T:xi)/Lx), 
namely it corresponds to having a single vortex threading 
the hole of the cylinder. Trivial momentum counting then 
tells us that this vortex carries momentum: Pvortcx — 

2'KVLy. 

Flux threading in the presence of boson interac- 
tions: 

In the presence of boson interactions, the phase ip at 
each site is no longer a c-number. Since the number and 
phase do not commute, interactions introduce phase fluc- 
tuations. Such fluctuations may permit the vortex to es- 
cape from the hole of the cylinder. Clearly, this is only 
possible if the initial and final state have the same crys- 
tal momentum quantum number. Since the vortex state 
carries Pvortox — "iiivLy, we expect the vortex will re- 
main trapped except at special values of v and Ly, where 
^vortex becomes a multiple of 27r. Thus, in general, even 
in the presence of boson interactions, threading flux $o 
introduces one vortex, carrying the above crystal momen- 
tum, into the hole of the cylinder. 
Flux threading and the Berry phase for a vortex: 

It is well-known that moving vortices in a station- 
ary Galilean invariant superfluid experience the so-called 
" Magnus force" , a force which acts transverse to the ve- 
locity of the vortex. A superfluid vortex thus behaves as 
a charged particle in a magnetic field, the Magnus force 



being analogous to the transverse Lorentz force. In a lat- 
tice system, this "magnetic field" seen by the vortex is 
encapsulated through vector potentials living on the links 
of the lattice, and the vortices pick up a Berry phase x 
(of the Aharonov-Bohm kind) on going around an ele- 
mentary plaquette of the lattice. We will now show that 
using the momentum balance argument fixes this Berry 
phase to be X = 27ri^. 

Consider a single vortex on an infinite plane. Since 
the vortex sees a "flux" x P^r plaquette of the lat- 
tice, the unit translation operators for the vortex satisfy: 
TxTy = TyTxexp{—ix)- Let \Kx,Y) represent the state 
with one vortex with x-crystal momentum Kx located at 
y — Y . When this vortex is translated by one unit along 
the -fy-direction, i.e. to y = Y+1, it is straightforward to 
show that the new state has x-crystal momentum given 
by Kx + X- For ^^ antivortex, the translation operators 
satisfy T^Ty — TyT^ exp(ix), and translating the antivor- 
ton along -|-y changes the crystal momentum to Kx — X- 

With this in mind, let us thread a vortex through the 
torus in the manner shown in Fig. ^ Start with a state 
with well-defined crystal momentum along x-direction, 
say zero. Create a vortex-antivortex pair on some pla- 
quette and make a superposition with zero net crystal 
momentum along x. Next drag them apart by translat- 
ing the vortex along +y using the translation operator Ty 
until they are L lattice spacings apart along the torus. 
This state then has additional crystal momentum xL- If 
we drag the pair all the way around the torus and annihi- 
late the vortex-antivortex pair, this would be equivalent 
to threading a vortex through the torus as in Fig.^b). 
The net momentum change is then x^y On the other 
hand, we have shown that Pf — Pi = 2'KvLy for thread- 
ing the vortex through the hole of the torus. This fixes 
X = 27r:^. 

We have confirmed this result by using the well-known 
duality mapping^- between bosons and vortices in 2 -|- 
1-dimensions. The dual theory treats vortices as point 
particles minimally coupled to a non-compact U(l) gauge 
field. In the dual theory, the flux of the gauge field on 
an elementary plaquette as seen by the vortex emerges 
naturally as x = 2'kv, i.e. the vortices see each boson 
as a source of 27r magnetic flux. The non-compactness 
of the gauge field, or the conservation of the magnetic 
fiux piercing the lattice, is a simple consequence of total 
boson number conservation. 

To summarize, vortices in a uniform superfiuid pick up 
a Berry phase xA/" on adiabatically going around a loop 
enclosing A/" plaquettes. The Berry phase per plaque- 
tte is completely determined by the particle density as 
X — 2^1/. Writing x — "^t^oim, which defines the "Mag- 
nus coefficient" q;m, leads to the Luttinger relation for 
superfiuids, namely 



I' = au + P 



(54) 



This relation follows from using the momentum balance 
arguments of Oshikawa, applied to a conventional super- 
fiuid. In this sense, the Berry phase relation above may 
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FIG. 4: Threading a vortex through the hole of a torus, (a) A 
vortex (arrow coming out) - antivortex (arrow going in) pair 
is created, and separated by the translation operator Ty. (b) 
When the vortex is taken all the way around the torus and 
then annihilated with the antivortex, the resulting state has 
one unit of circulation about the hole of the torus as shown. 



be viewed as the analogue of the Luttinger relation for 
Fermi liquids. 

There is a concern which we have not addressed so far 
— the vortex may have a modified density near its core, 
and this in turn could modify the Berry phase accumu- 
lated by the vortex when it is adiabatically taken around 
a loop. However, this term does not change with the 
area of the loop, so we can still use the above result to 
deduce a precise difference of the Berry phase between 
two loops enclosing different areas. We need to define 
2TiaM = A^/AA/", where A$ is the difference in Berry 
phase between two loops which differ in area by ZV^/' pla- 
quettes. Defined in this manner, 2'KaM obeys the precise 
relation in 1)54(1 . Another caveat is the very definition of 
adiabaticity in the presence of gapless superfluid phonon 
excitations — to be completely rigorous, we need to work 
with a finite-sized system such that the phonons have a 
nonzero gap, and demand that the vortex motion be adi- 
abatic with respect to this energy scale. 



X. Z2 FRACTIONALIZED SUPERFLUID ST* 



bosons. This supports two gapped excitations: charge 
(5/2 chargons that also carry Z2 gauge charge and Ising 
vortices (visons). On doping this insulator, the addi- 
tional charge-Q bosons can deconfine into pairs of char- 
gons since in the fractionalized phase. Bose condensing 
the doped chargons would destroy deconfinement of the 
Z2 gauge field (by the Anderson-Higgs mechanism, since 
the condensate carries Z2 gauge charge) and lead to a 
conventional superfluid. 

The other possibility is that doped chargons pair and 
Bose condense resulting in a superfluid phase. Since the 
condensate is Z2 gauge neutral, deconfinement is pre- 
served and this exotic superfluid phase is called SJ^* . ft 
supports elementary hc/Q vortices, and the visons still 
survive in the superfluid. There is however another ex- 
citation, analogous to a Bogoliubov quasiparticle of a 
superconductor, present in the system — this gapped 
quasiparticle^ is the "ison" . It may be viewed as a de- 
scendant of the chargon in the insulator, whose electric 
charge has been screened by the Bose condensate of pairs 
so that it only carries a Z2 Ising chargeSS. In addition to 
these gapped excitations, there is of course the gapless 
superfluid phonon in an electrically neutral system. The 
relative statistics of the gapped excitations are as follows: 
the wavefunction changes sign if an ison is adiabatically 
taken around the vison or the vorton. 

We have seen how Z2 fractionalized Bose Mott insula- 
tors with full lattice translation symmetry fall into two 
classes, depending on whether the boson flUing is an in- 
teger (I|i,e„) or half odd-integer (T*^^ as described in 
Section El Similarly, fractionalized Fermi liquids !FC* 
also come in two varieties as shown before, with different 
relations between fermion filling and Fermi volumes. It 
is not surprising therefore that we will find below two 
kinds of Z2 fractionalized superfluids, ST* . Using mo- 
mentum counting arguments as done earlier, we will also 
see how the distinction between the two types of ST* 
phases, namely STly^„ and ST*^^^^, is reflected in differ- 
ent dynamics for the vorton in these two cases. While in 
the case of translationally symmetric Z2 insulators T* , a 
knowledge of the flUing alone was enough to determine 
the odd/even nature of the phase, this is no longer true in 
the case of Z2 fractionalized superfluids, ST* (or TC*), 
where additional information regarding the odd/even na- 
ture of the phase is needed. 



We now discuss an exotic variant of the superfluid, 
ST*, and see how the relation H54|) is modified in this 
phase. The ST* phase is a Z2 fractionalized superfiuid 
which was first discussed by Senthil and Fisher.^ It sup- 
ports three distinct gapped excitations: (i) an elemen- 
tary hc/Q vortex (called the "vorton" in Ref. W), (ii) Z2 
gauge fiux (the "vison" ) , and (iii) an electromagnetically 
neutral particle carrying Z2 gauge charge (the "ison"). 
There are various ways in which the ST* phase may be 
realized, we shall briefiy outline one of them. 

Let us start with a fractionalized insulator X* which 
can be realized at integer or half-odd integer density of 



A. Effective Hamiltonian for ST* 

A simple description of the ST* phase may be ob- 
tained by using a Hamiltonian which describes chargons 
minimally coupled to a Z2 gauge field. This takes the 
form 



with 






H{ST*) = f+V + Hg (55) 

QAj/2) 



ij COS((pi — (pj -r ^j-i-i'j I 
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- tB^cos{2(t)i-2(l)j + QA,j) 

V = MntN 



(u) 



(56) 

(57) 



and physical states of the theory satisfy the constraint 
in Eqn. 1211 Here exp(— z0i) creates a chargon carrying 
charge Q/2 and Z2 gauge charge at site i, rii is the char- 
gon number, and the terms in T represent the chargon 
kinetic energy. Singie-chargons hop with an amphtude 
tf), and are coupled minimally to the Z2 gauge field and 
for simplicity of discussion, we have included explicitly a 
chargon-pair hopping term with amplitude ts- Clearly 
a chargon-pair created by exp(— 2i0i) has no net Ising 
charge and does not couple to the Z2 gauge field. T^nt [n] 
is an interaction term involving chargon densities which 
we do not spell out here. The exotic superfluid phase 
SJ^* requires being in the deconfined phase of the gauge 
theory (which is guaranteed by a large K/h) and with 
chargon pairs condensed (which can be achieved with 
large chargon pair hopping ts , while single chargon hop- 
ping remains small). 

Let us now write down the effective Hamiltonian in the 
5jF* phase. Condensation of chargon pairs implies that 
we can replace the operator exp(— 2i0i) by a c-number. 
Then, the magnitude of the chargon creation operator is 
determined, but its sign can fluctuate which gives rise to 
the ison field - i.e. we can write exp(— i^^) c>c If, where 
If is a Pauli matrix with eigenvalues ±1. Similarly, since 
chargon-pairs are condensed, the parity of the chargon 
number operator chargon number n^ must be changed 
by the ison creation operator If , hence we identify rii w 
(1 -|- -ff )/2; again /^ is a Pauli matrix and (1 -I- /^)/2, 
with eigenvalues {0,1}, counts the number of unpaired 
Ising charged particles. 

In new variables, the Hamiltonian (with Aij = 0) re- 
duces to 



Hic 



Tied + K-cd + Hg + -^condensate (58) 



-t: 



' '^ IfP.CJ- 



Tied 



E 



(59) 



where we have introduced a "chemical potential" g for 
the isons, and -ffcondonsato describes the dynamics of the 
condensate. Physical states of this theory need to satisfy 
the constraint: 



n 



-If 



(60) 



^(^) 



This is simply an Ising model coupled to an Ising gauge 
field - the ST* phase is realised when both the isons 



(excitations of the Ising model) and visons (excitations of 
the Ising gauge theory) are gapped. This will occur when 
K/h and l^/i^l are large, when the gauge theory is in the 
deconfined phase and the Ising model is 'disordered'. Let 
us now briefly consider the special limiting cases where 
\g\ ^ cx) in order to expose the underlying reason for the 
two kinds of Z2 fractionalized S!F* phases. Clearly, if 
g — > ±00, we would have /f = ±1 corrresponding to the 
ison number (1 + Ti)/2 = 0, 1 respectively. The physical 
states of the gauge theory then satisfy the constraint: 
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=Foo 



: n 

i=nn(i) 



a^=±l 



(61) 



These constraints on the gauge theory, as we know 
from the discussion on insulators, correspond to "even" 
and "odd" Ising gauge theories respectively, which corre- 
spond to having zero or one Ising charged particle (ison) 
fixed at each site. Thus, for g — > -f-cxD and for g — > —00 
we will obtain two distinct superfluid phases, which we 
label SJ-ly^^ and SJ-*^^^ respectively, which persist to 
finite values of g as well . These are separated by an in- 
termediate phase where \g\ ^ t^ where the /^ Ising field 
orders; this is the conventional superfluid phase. Below, 
we will see how these ST* phases can be distinguished 
from each other. 



B. Flux threading in ST* 

Trivial momentum counting tells us that threading flux 
$0 = hc/Q into the cylinder on which the system lives 
changes the crystal momentum by 2'KvLy, where v is 
the filling and Ly is the number of rows of the cylin- 
der. Where is this momentum soaked up in the ST* 
phase? We will show below that flux threading intro- 
duces both a vison and a vorton into the hole of the 
cylinder. The crystal momentum is then divided up be- 
tween these two excitations, in a way that depends on 
whether we are dealing with ST*^^^^ or ST*^^^. This will 
be argued below in two stages. First, we consider freez- 
ing the Ising gauge field dynamics by setting the vison 
hopping to zero (i = 0). There it can easily be argued 
that $0 flux threading leads to both a vison and a vorton. 
Then, using our earlier knowledge of vison momenta in 
the even/odd Z2 gauge theories, and the total momen- 
tum imparted to the system, we can read off the crystal 
momentum carried by the vorton. Finally, we reinstate 
the vison hopping [t > 0) and use continuity to argue 
that this does not affect the momentum assignments. 

Threading flux (f>o introduces a vison and vorton: 

Consider at first the limit h ~ and ts — 00 identically, 
so that before introducing flux we can everywhere set 
cr? — 1 (as a reference state which we can then project 
into the subspace of physical states) and exp(z2(/)i) = 1. 
Let us adiabatically thread flux hc/Q in the —y di- 
rection for the above system on a cylinder as in Fig. ^ 
such that starting from the initial eigenstate |^(0)) in the 
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absence of flux, the final state reached is |^(T)). The fi- 
nal state can of course be written as |^(T)) = Ut\'^{0)) 

with Ut = Ttexp (— j/q Ha{SJ-* ,t)dt\ where Tt is the 

time-ordering operator. We can go to the Aij — gauge 
by making a unitary transformation Ha {ST* , T) — > 
UgHa{SJ-* ,T)Uq^ = Hq{ST*) (corresponding to zero 
vector potential). Here Ug = U^Ua-, with 



Ua 



U„ 






(62) 



ijGcut 



P alj where the loop C is taken around the 



and 'cut' refers to the vertical column of links for which 
Xi — Lx,Xj — 1 (shown in Fig. |2Ib)). Thus, the fi- 
nal state in the Aij = gauge is |^/) = U(I,UMt\'^ i) ■ 
Since the system is a superfluid initially in the state 
I exp(2i(/)i) — 1), which is unaffected by the time devel- 
opment operator Ut since we are in the ts — oo limit, 
the final state after acting with U^ has | exp{2i(f>i) — 
eicp{i2TTXi) / L^) , namely it corresponds to having a single 
vorton threading the hole of the cylinder. 

At the same time, following arguments similar to the 
insulator X*, acting with U^ introduces a vison in the 
hole, the vison number being defined by (1 — Wc)/^ with 

Wc = U, 

cylinder (see Fig.OIb)) 

Thus, for h = 0,tB = oo, threading a 27r flux adds a 
vison and a vorton into the hole of the cylinder. Thus, 
when the effective description of the gauge fields is an 
even (odd) Ising gauge theory (which are connected to 
the g — > — oo (+oo) limits respectively as we have seen 
above), momentum counting arguments already showed 
us that the vison carries momentum zero (for the even 
gauge theory) or nLy (for the odd gauge theory). The 
remaining momentum must clearly be carried by the vor- 
ton! 

Thus, in the case of even Z2 gauge theories the vison 
carries zero crystal momentum, and we deduce that the 
vorton through the hole of the cylinder carries crystal mo- 
mentum 2TTh'Ly, and this is the ST*^^^^^ case. For the case 
of odd Z2 gauge theories , the vison carries crystal mo- 
mentum T^Ly and the vorton carries crystal momentum 
2'k{v— l/2)Ly, and this is the ST*^^^ case. To summarize, 
momentum balance arguments suggest 



In this case the loop product Wc no longer commutes 
with the Hamiltonian, we cannot use its eigenvalues to 
label the eigenstates. Thus, in a finite system as argued 
previously for the case of inulators, the vison can tunnel 
out in the case of even gauge theories for any cylinder 
dimension, or in the odd gauge theories, if the cylinder 
has an even number of rows {Ly). However for the case 
of odd gauge theories on a cylinder with odd number of 
rows {Ly), the vison carried momentum tt, and hence vi- 
son tunneling is blocked even in finite systems. Similarly, 
when there is dynamics to the condensate (is ^ 00), the 
vorticity is also not necessarily conserved in a finite sys- 
tem, namely / V2(/) around the cylinder is not a constant 
(not a classical variable). More precisely, this statement 
can be made in terms of an order of limits; if the time 
scale for flux threading is taken to infinity before the 
thermodynamic limit is taken, then the system can re- 
main in the zero vorticity state at the end of the flux 
threading. This occurs if the two states, namely the 
state with no vison and no vorton and the state with 1- 
vison and 1-vorton, each carry zero crystal momentum, 
and can then mix to give eigenstates of the Hamiltonian. 
This happens if 2'KvLy = O(mod 27r). Otherwise, even in 
the presence of h and condensate dynamics, the vorton 
acquires a nonzero crystal momentum and therefore its 
tunneling is blocked even in a finite size system (in the 
sense described above). Thus, in all cases where a non 
trivial crystal momentum is imparted to the system from 
flux threading, this is accounted for by the presence of a 
vorton and/or a vison in the flnal state which carries the 
appropriate momentum. 

The two superfluids SJ-*^^^^ , ST*^^^ can then be distin- 
guished depending on the momentum carried separately 
by the vison and the vorton though the total momentum 
carried by these excitations is the same. We shall later 
see that this may be reflected in the measured Hall effect 
in the vorton liquid phase in these systems through the 
Berry phase induced Magnus force on the vorton. In the 
next subsection, we shall directly verify the momentum 
assigments for the vorton by going to a dual theory where 
vortices are represented as particles. 



C. Berry phase for vorton and consistency with 
momentum counting 



rteven 
vorton 

nodd 



2TriyLy[niod 27r] 



vorton = 27r(zy - -)ij^[mod 27r] 



(63) 
(64) 



The dual theory for the ST* phase, where the chargons 
are traded for vortex variables, is derived in Appendix A. 
The action takes the form 



We now show that these momentum assignments are not 
affected on reinstating the vison hopping {h ^ 0) and the 
condensate dynamics (t^ ^ cxd). 

Flux threading with dynamical gauge fields: Here, 

starting from the case with h = 0,tB = cxd, let us ask 
what happens if we turn on a nonzero h and finite t^, 
giving dynamics to the gauge field and to the condensate. 
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+ c^ J2 iJi.y+^Uj2{J^~nr. (65) 

z,/i— 1,2 i 

Here exp(— i0) creates a vorton, and the first term in 
Sc represents vorton hopping. The vortons appear as 
charged particles minimally coupled to the gauge fields 
Af_i and Gfj,. The flux of the (noncompact) U{1) gauge 
field A^ and the gauge field a^ are tied to the charge 



current (J7^) and the ison current (7^) respectively: 



^7-^/2 









(66) 
(67) 



The isons have a chemical potential g, and a Berry phase 
term associated with the fact that they couple to the 
Z2 gauge field tXy . The remaining terms represent local 
interactions between the charge density /currents. 

On the spatial links (/x = 1, 2), (/^) = and (J^) = 0. 
On the temporal links, for the chargon density we have 
(J'q) = 2^, where v is the average charge density in units 
of Q (equivalently, i' is the density of chargon-pairs) . 
For the ison density, we have two possibilities: (i) for 
g — > — cx), we have /q = and there are no isons in the 
ground state; (ii) for g —^ +cx3, there is one ison nailed 
down to each lattice site, Tq ~ I. For g -^ —00 (or 
5J-"*„g„), when the vortons go anticlockwise around an 
elementary plaquette of the square lattice, they see only 
the flux produced by the vector potential 2nA, and the 
wavefunction acquires a factor exp^inJ'Q). On average, 
the phase picked up is thus 27rz^. This is identical to the 
Berry phase picked up by a vortex in a conventional su- 
perfluid with boson density v. For g —^ +00 (or SJ^'^^^), 
the vorton sees an additional flux produced by one ison 
charge Iq — I nailed down to each site, and the total flux 
seen by the vorton is thus 27r(z/— 1/2) — this deviates by 
TT from the Berry phase for 5J^*^g„ and the conventional 
superfluid. 

Applying the picture of vortex threading presented for 
conventional superfluids to this case of vorton threading, 
it is clear that this Berry phase is consistent with the mo- 
mentum counting argument. Namely, the vorton thread- 
ing suggests that x^y = i^vorton- On the other hand, 
momentum balance tells us Pf — Pi = iirvLy for thread- 
ing a vison and a vorton through the hole of the torus. 
Since we know that the vison carries crystal momentum 
in an even gauge theory (as in 5J-"*^g„) or 
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^) _ 



vison 
(odd) 



-Wison — ^^y 



in an odd gauge theory (as in ST'^^^). 
This fixes Pvorton in the two cases (in agreement with 
Eqns. IMEl and thus x*™^"^ = 2Tn^ (in SJ^l^^J, and 
^{odd) ^ 27r(iy- 1/2) (in STl^^). This is consistent with 
the result derived from the dual theory above. Defin- 
ing as before the "Magnus coefficient" aj\/ = x/^tt, we 
obtain 



„ even 



„ odd 



p 



{^ 



\^ 



p, 



(68) 
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FIG. 5: Schematic figure showing the filling dependence 
[v = (1) is the empty (full) band) of the Hall conductivity 
in a conventional Fermi liquid or TCtve.n (solid line), for elec- 
trons within a simple Drude-like picture, contrasted with the 
behavior expected in the TC^^^ phase (dashed line). Similar 
results are expected for the conventional superfiuid or SJ-'^ven 
phase compared to the SJ-^^^ phase in a vortex liquid phase. 



where p is an arbitrary integer. Thus the momentum 
counting provides a prescription to fix the Berry phase 
for the vorton, and allows us to distinguish the "odd" 
and "even" exotic superfluid phases. 



XI. CONCLUSIONS 

Extending a non-perturbative argument, made by Os- 
hikawa for the Fermi liquid, we have constructed ana- 
logues of Luttinger's theorem for systems other than the 
conventional Fermi liquid in dimensions D > 2. This 
has allowed us to derive constraints which must be sat- 
isfied by quantum phases of matter on a lattice, such as 
superfluids and the more exotic Z2 fractionalized phases 
which are topologically ordered. We have discussed ways 
in which these constraints may be useful in identifying 
fractionalized phases in numerical experiments. 

A recurring theme has been the important distinction 
between 'even' and 'odd' deconfined Ising gauge theories, 
which correspond to states that are most naturally asso- 
ciated with integer and half integer filled systems respec- 
tively. For exotic insulators, the even or odd character 
of the phase is completely determined by the filling in 
this manner. For exotic Fermi liquids and superfluids, a 
knowledge of the fllling by itself is insufficient to deter- 
mine the odd or even nature of the emergent Z2 gauge 
field — precise violations of the Luttinger relation or its 
analogue for these systems provides a way to distinguish 
them from each other. 

Within a simple Drude-like picture, one associates the 
size of the Fermi surface to the sign of the Hall conduc- 
tivity — a Fermi surface corresponding to a few electrons 
exhibits an electron-like response, while a Fermi surface 



23 



of a nearly filled band would show a hole-like response. A 
Fermi surface which violates the conventional Luttinger 
theorem may thus be reflected in an anomalous sign of 
the Hall conductivity as depicted schematically in Fig. O 
A similar change in the sign of the Hall effect in a vortex 
liquid phase is expected for odd fractionalizcd supcrfluids 
relative to conventional superfluids, due to a shift of the 
Berry phase by tt at a given density of bosons. Clearly, 
the sign of the Hall effect is not universal and in real 
systems is affected by band structure and interactions. 
Thus an anomalous sign of the Hall response is sugges- 
tive but is not a rigorous diagnostic. Experimental tools 
such as angle resolved photoemission spectroscopy which 
measure the Fermi surface can more directly detect viola- 
tions of the conventional Luttinger theorem expected in 
odd Fermi liquids and thus serve to identify such systems. 
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APPENDIX A: DUALITY AND VORTONS IN ST* 

Let us begin with the path integral for the partition 
function of chargons coupled to a Z2 gauge field, Z = 
/^</'S{„}S{o^}exp(-S'), with S = Sgauge + Sc- The 
gauge field and chargon actions are given by 



-sKs Y[ CTij - eKr W^ (Tij 

{ijra i 

'.'^ni{4>i - (j)i+r + -^[l - (^iS+r]) (Al) 



Here, (ij) denotes nearest neighbor sites in space, n, 
denote the chargon number and phase. The chargons 
hop with an amplitude if, and have a local repulsion of 
strength U . Un plays the role of the chargon chemical 
potential. The Z2 gauge fields are denoted by (7^ — 
±1 and Dj,, D.^ denote elementary spatial/spatiotemporal 
plaquettes on the cubic space-time lattice with respective 
gauge field couphngs Ks,Kr- Finally, e is the Trotter 
discretization along the imaginary time direction. 

We can rewrite the chargon hopping term in the par- 
tition function as 



Lii 



,]) 



(A2) 



where Lij = —Lji is an integer-valued field. For large 
a, with a — ln(2/etb), this reduces to the original char- 
gon hopping term. For general a, this modified form 



allows terms such as cos{2(j)i — 2(pj) which correspond to 
chargon-pair hopping. We therefore do not need to keep 
an explicit pair-hopping term ts unlike in our discussion 
in section Ixl 

Integrating out the phase field <j> leads to a constraint 



y ^ Lij + {rii - rii+r) = 



(A3) 



which is just the discrete version of the continuity equa- 
tion, Lij representing the chargon current on bond (z — > 
i). Writing the number n^ 
in the form 



Lii+Ti we can recast this 



/j=0,l,2 



i.i+ajj. 



(A4) 



where xo = r, xi = x, X2 — y- This sets the divergence 
of the 3-current to zero. Below, the sum over /i will be 
understood to run over 0, 1, 2 unless stated, we will also 
use the notation L' = Lii^^,^ ■ 

We go to dual vortex variables in 2 -I- 1 dimensions 
in the standard manner^^, the only difference is in the 
presence of Z2 gauge fields in the action but we do not 
dualize these. The constraint is solved by equating the 
conserved current to the curl of a dual vector such that 
its divergence is automatically zero. We decompose the 
chargon current into two parts, the current of pairs J 
(an even integer) and the current of unpaired particles / 
{— 0, 1). Note that / is only conserved modulo-2 — two 
unpaired particles can combine to form a pair which is 
accommodated by increasing J by one unit, and / thus 
is the current of particles carrying only a Z2 charge. The 
constraint is thus solved by choosing 



11 = (e^'-^a.aDimod 2], 



(A5) 
(A6) 



where A (an integer) and a (= 0, 1) arc fields on links 
of the dual space-time lattice, and the right hand sides 
above are just the lattice curls on the dual lattice, taken 
around the original link (i, i -I- x^). 

The chargon action then takes the form 

Sc^ a Y. {Jl + Ilf + eU{Jl, + r,~nf (A7) 
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where we have now included a chemical potential g for 
the Z2 charges whose density is Iq = e^^^dvaxlmod 2]. 

We can convert the sum over A to an inte- 
gral by softening the constraint by introducing terms 
£it,^j cos(27rA^) in the action (this can be formally 
accomplished by using Poisson resummation), which 
prefers A^ to be an integer. Everywhere else in the ac- 
tion, only the transverse part of A plays a role (since 



only its lattice-curl appears). Extracting the longitudi- 
nal part of 2nAij as 9i — 9j , we identify the dual vorton 
creation operator exp{—i9i). The vortons are seen to be 
minimally coupled to the transverse part of the A, which 
we denote A, exactly as a charged particle coupled to a 
U{1) gauge field. Thus, in the softened theory, 



27rA! 



(^i+x„ 



+ 27r^j,. 



(AlO) 



Making this substitution, and absorbing a by shifting 

■K^-K- «^/2 find 



E (^;)'+E('^o-")' 



(All) 
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'l,fj, 



t, E cos(0. - 0,+,^ - 2tiA'^ ~ ^a^) (A13) 



with the total current J' = Jl + 11 = 2e„^A9^-4\. This 



is the result used in Section IXIC 
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